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SECTIOH I 


INTRODUCTION 


NASA Dryden Flight Research Center has the responsibility for flight 
test of advanced remotely piloted research vehicles (RPRV) to explore 
highly maneuverable aircraft technology, and to test advanced structural 
concepts, and related aeronautical technologies which can yield impor- 
tant research results with significant cost benefits. These RPRV's, 
such as HIMAT (Ref. 1) and DAST, are launched from a B— 52 and controlled 
remotely from a ground pilot station. Fig. 1, at DFRF using a telemetry 
uplink for issuing control commands and a downlink for sending sensor 
data. The RPRV is tracked by a ground radar station capable of tracking 
the vehicle with less than 12 yards error. The radar processing facil- 
ity is capable of generating guidance information over the telemetry 
uplink including synthetic ILS localizer and glide slope data for three 
runways on the lake bed. 

At the present time landings on one of the three runways is accom- 
plished by visual reference to a forward looking TV image which is data 
linked to the remote cockpit CRT display. Prior to landing flare the 
pilot has synthetic glide slope and localizer errors displayed on the 
remote cockpit attitude director indicator (ADI) display. Using the 
ground station control stick and rudder pedals the pilot flys the RPRV 
to acquire the synthetic beam and stabilize the vehicle path along the 
approach path by minimizing the director indicator errors. For the 
flare and touchdown the pilot uses the TV image for landing. The pilot 
lacks altitude data and a yaw angle due to wind crab looks like a runway 
offset. Consequently, the last 500 ft of landing is confusing. The 
pilot work load is high, and difficult to accomplish, but the pilots do 
an acceptable job after performing several practice landings in a PA-30 
aircraft with a safety pilot prior to an actual RPRV flight. 

An automatic flare mode is available for the longitudinal axis as a 
part of the HIMAT backup control system. This mode can be activated 
when below 5000 ft radar altitude by selecting the LAND mode. The 
vehicle then executes a preprogrammed flare and airspeed reduction 
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maneuver referenced to radar altitude but again the T.V. display to the 
pilot is confusing* 

In case of a failure in the uplink the RPRV is commanded into a 
level, constant bank angle orbit which continues until communication is 
restored or the vehicle crashes* At the present time there is no means 
to command the vehicle to automatically proceed to a safe unpopulated 
area in case of total uplink loss* 

NASA DFRF is interested in developing a system which shall make man- 
ual landings easier and more consistent and to provide a simple autono- 
mous navigation system aboard the RPRV to have it proceed to a safe area 
for orbiting after loss of uplink. The provision of an autoland capa- 
bility supervised by the remotely located pilot would be beneficial in 
reducing pilot workload and making the landings more consistent. 

The primary purpose of this study is to provide the preliminary 
design of an upgraded automatic approach and landing control system and 
flight director display to Improve landing performance and reduce pilot 
workload* A secondary purpose is to determine the feasibility of an 
onboard autonomous navigation, orbit, and landing capability for safe 
vehicle recovery in the event of loss of telemetry uplink communication 
with the vehicle. 

The following secti^ons . describe the current RPRV approach and 
landing method, the proposed automatic and manual approach and autoland 
system, and an autonomous navigation, orbit, and landing system concept 
which is based upon existing operational technology* 

The intent of the study is to upgrade the DFRF facility rather than 
a specific vehicle performance. At the outset of the study the HiMAT 
RPRV was specified as the example vehicle. However, it rapidly became 
apparent that the desired flight control and director display algorithms 
must be vehicle specific and the HiMAT vehicle would no longer be opera- 
tional by the time any of the control and display concepts could be 
implemented. Since the PA-30 aircraft has been used in the past as a 
developmental test bed and is likely to continue in this role in the 
future, the various analysis herein have been more or less focused on 
the PA-30. 
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SECTION II 


OJRRENT RPRV APPROACH AND LANDING 

A. MANUAL APPROACH AND LANDING 
1 • Normal Procedures 

The approach and landing is flown by the ground based pilot viewing 
the head down instrument panel or head up TV display and operating the 
cockpit controls (stick and pedals). The pilot is guided to intercept 
an extension of the runway centerline via verbal guidance commands of 
test personnel who monitor the horizontal position of the vehicle on an 
X-Y plot board containing a map of the area. The vehicle is then 
descended, at a rate based upon the specific test conditions, to inter- 
cept a computer generated "synthetic ILS" glide slope of 2.82 deg emi- 
nating from a runway aim point one nautical mile from the runway thres- 
hold (see Fig. 2). 

Measures representative of deviation from the "synthetic ILS" beam 
are calculated in the 620~f-100 computer based upon FPS-16 tracking 
radar data (see Fig. 1). These deviations are modified in V-73 Compu- 
ter A to drive the flight director needles of the Attitude Director 
Indicator (ADI) on the pilot's display panel. Thus, the pilot can "fly" 
the "synthetic ILS" beam until the runway comes into view of the vehicle 
mounted TV camera and is displayed on the head-up TV monitor. The final 
approach and flare is then flown by means of the TV view with an engi- 
neer calling out altitude and airspeed. 

Stick and rudder pedal deflections by the pilot command pitch, roll, 
and yaw rate responses via the control algorithms mechanized in the 
ground computer. The computer outputs surface actuator position com- 
mands which are then uplinked to the RPRV. 
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Figure 2. Typical HIMAT RPRV Approach Paths 


























2m ADI Synthetic ILS laformatlon 

Documentation was not found for the algorithms used to compute the 
synthetic ILS beam and vehicle displacement error signals used in the 
ADI flight director mode. However, error could be determined from radar 
tracking measurements and trignometric relationships such as depicted in 
Fig. 3. The FPS-16 tracking radar provides measurements of the RPRV 
range, R, azimuth, A^, and elevation, E^, relative to the tracking sta- 
tion (or some other bench mark). The RPRV location relative to the syn- 
thetic beam is computed to produce glide slope and localizer position 
error measures, h^ and y^, respectively to drive the ADI needles. 

Based upon the relationships shown in Fig. 3 the error calculations 
could be: 

y^ = Y - R cos S 3 
h^ ^ h h^ 

^ (R sin £3 4* Z) - L tan y 

where 

Y = perpe'ndicular distance from runway center- 
line to radar 

Z = height of radar above like bed 

L = R cos Cq (sin - cos tan af) 

R = measured vehicle range from radar 

= measured vehicle elevation above radar 

£3 = measured vehicle azimuth from perpendicu- 

lar to runway 

= angle between perpendicular to runway and 
azimuth to aim point 
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Figure 3. Pseudo or Synthetic ILS Beam 


The FPS-16 measurement accuracy is (Ref* 2): 

Range: i30 ft 

Azimuth: ±0*1 mil 

Elevation: ±0.1 mil 

The radar is located on a hilltop approximately 2.5 nm west of the 
intersections of runways 15/33 and 5/23 (see Fig. 4). From this vantage 
point it provides an unobstructed view of these runways and runway 7/25, 
the runways designated for RPRV landing (Ref. 3). Unfortunately, the 
orientation between the runway centerline and radar location is such 
that vehicle lateral (localizer) error would probably be determined from 
the range measurement on runways 15/33, 5/25, and 5, and from azimuth 
only on runway 23. Although this has some influence on error measure- 
ment accuracies the impact is small considering the width of these run- 
ways and the surrounding areas of lakebed. The principal effect would 
be on the algorithm to compute lateral error. 

3. TV Display 

The forward looking TV cameras provide a restricted, two dimensional 
black and white view of the area directly ahead of the vehicle. The 
camera axis is depressed with respect to the aircraft longitudinal axis 
to assure a view of the^ runway during the flare. Thus, in steep des- 
cents the horizon is not visible and roll attitude information is lack- 
ing. The runway does not come into view until late in the approach and 
then it is difficult for the pilot to extract the necessary information 
as to runway alignment, flight path, altitude, etc. As an example, 
three different restricted views of the runway from a banked aircraft 
are shown in Fig. 5 (from Ref. 4) along with an unrestricted view. All 
views are the same runway and horizon sketch but only in the unrestric- 
ted view can it be determined that the aircraft is aligned with the left 
side of the runway. 

The monitor, an 11 in. (diagonal) black and white display, is 
located about 20 in. from the pilots eye. The camera has a 50 ram lens. 
This results in a magnification factor of 1.8 (Ref. 4). 
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Figure 4. Runways Designated for RPRV Landings 
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During the final approach and flare phase the pilot is reluctant to 
look away from this restricted view long enough to obtain altitude or 
airspeed data from the instrument panel* Thus, the necessity for an 
engineer to call out this information* 

4* Control System 

The remote vehicle control mechanization described here is the sys- 
tem implemented for the PA-30 aircraft (Refs. 5 and 6) to validate the 
HIMAT backup control system (BCS) autoflare concept prior to the first 
HIMAT flight* This PA-30 mechanization subsequently has been employed 
as a ground pilot trainer/ref resher system prior to each RPRV flight. 
In this application the RPRV flight control system is mechanized on the 
ground computer and commands to the various control actuators are 
uplinked to the PA-30 servo amplifier* The vehicle motion parameters 
downlinked to the computers are listed in Table 1 (from Ref* 5)* 

The PA-30 has an onboard pilot to fly it via its own control system 
to the desired area for initiating the ground controlled approach. Con- 
trol is handed off to the ground pilot in what is termed the computer 
direct (CD) mode for pitch, roll, and yaw* See Figs* 6, 7, and 8 taken 
from Ref. 5* The ground pilot then selects the pitch, roll, and yaw 
rate damping (RD) modes as the primary means of maneuvering control. 
This mechanization incorporates both rate and attitude feedbacks in the 
pitch and roll axes and are a form of rate command, attitude hold con- 
trol* The latter modes had been used previously in the PA-30 aircraft 
and gave improved handling over the basic control system* 

When approach line-up requirements have been met the ground pilot 
can select the autoflare mode for completion of the approach and 
landing. 

B. HIMAT BACKUP CONTROL SYSTEM ADTOFLARE MODE 

The ”Autoflare** mode is pilot selectable at any altitude below the 
radar altimeter upper limit of 5000 ft* The longitudinal pitch axis 
incorporates a preprogrammed two-step altitude rate descent followed by 
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TABLE 1. PA-30 DOWNLINKED VARIABLES 



MNEMONIC 

DESCRIPTION 

SCALING 

1 

P 

Roll Rate 

±40 deg/sec 

2 

Q 

Pitch Rate 

±25 deg/sec 

3 

R 

Yaw Rate 

±40 deg/sec 

4 

Theta 

Pitch Attitude 

±60 deg 

5 

Phi 

Roll Attitude 

±90 deg 

6 

HRAD 

Radar Altimeter 

0-5000 ft 

7 

DWl 

Digital Word #1 

— 

8 . 

DW2 

Digital Word #2 

— 

9 

VKIAS 

Indicated Airspeed 

0-200 kts 

10 

SQSWM 

L/H Gear Squat Switch 

On/Off 

11 

i DAP* 

1 

L/H Lateral Wheel Position 

±10.2 in. 

12 

1 

DRP* 

L/H Pedal Position 

±5.24 in. 

13 

DEP* 

L/H Pitch Wheel Position 

±13.16 in. 

14 

A2 

Normal Accelerometer 

+2 to -5 g 

15 

AY 

j Lateral Accelerometer 


16 

HOOT 

j 

IVSI 

±2000 ft/min 


During RPRV mode of operation the ground pilot control positions 
will be used instead of the variables marked with the asterisk. 


TR- 1180-1 


12 






DETRM 



TR-1 180-1 


13 
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0.5 

(deg /sec) 
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Figure 8. PA-30 Yaw CD/RD 






an exponential flare to landing* The complete approach and flare is 
scheduled via onboard radar measurement of altitude above ground level 
but provision is made for the pilot to insert discrete climb/dive rate 
increment corrections if necessary. The longitudinal axis also contains 
an autothrottle for airspeed control — again scheduled via radar 
altitude. 

The roll axis '*autof lare” mode is a zero yaw rate hold system with 
pilot coramandable discrete left/right step turn rate commands. The yaw 
axis consists of a yaw rate damper with lateral acceleration turn 
coordination. 

!• Longitudinal Mechanization and Performance 

PA-30 pitch axis and autothrottle mechanization block diagrams are 
presented in Figs. 9 and 10, respectively. In the pitch axis, the path 
sink rate control is implemented through the and hg^Q feedback 
loops. The inner loop normal acceleration and pitch rate feedbacks pro- 
vide . the proper equalization and vehicle damping for the outer path loop 
control. A stability analysis of the PA-30 pitch feedbacks including 
sensor and actuator dynamics is contained in Appendix A. 

To provide a reference for the improved autoland system proposed 
later, an idealized PA-30 system performance analysis will be made in 
which sensor dynamics are assumed to be perfect and the sink rate and 
airspeed schedules are approximated as shown in Fig. 11. These com- 
manded schedules divide the autoflare program into three phases which 
are summarized in Table 2. For the present assessment, only the decel- 
eration and flare regions are considered. 

Considering the altitude control channel as shown in Fig. 9 the sink 

« 

rate error, h^, is 


= “*^BARO + he 


( 1 ) 
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Figure 9. PA-30 Pitch Autoland Mechanization 
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- 16.6 


High-Pass Log 



Figure 10, PA-30 Autothrottle Mechanization 
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TABLE 2. THREE PHASES OF PEIESENT BCS AUTOFLARE 
AS MECHANIZED ON PA-30 




constant speed descent region 


hR > 950 ft 
h = -1000 ft/min 

= 110 kts (airspeed) 


(B) — deceleration region 
.150 ft < hR < 950 ft 
h = -500 ft/min 

decreases linearly with from 110 ft to 80 kt 

@ - flare regi on 

0 < Hr < 150 ft 

|hj decreases linearly with from -500 ft/min to -45 ft/min 
Va = 80 kt 
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-1000 


a> 

o 

^ c 



Radar Attitude, hp{ft) 

b) 


Figure 11* Simplified Sink Rate and Airspeed Schedules 
If the system works perfectly, the sink rate error is near zero, i*e 
he = 0 

Thus, the sink rate closely follows the h schedule 


hfiARO = = -500 ft/min = -8.333 fps 


Since the sink rate is constant in the deceleration region, 
altitude is 
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= f H 

•'O 


dt + h 


D 


~ hj^ + ht 


= 950 - 8*333t ft (3) 

where t = 0 @ = 950 ft. Thus, over the deceleration region from 

950 ft to 150 ft, the altitude change. Ah, is 


Ah = 150 - 950 = -8.333 to 


(4) 


from which 


to = 96.0 sec 

at start of flare (independent of wind). 

The throttle loop may be considered in a similar fashion resulting 
in 


Va 



= 125.7 -h 0.0634 h^ fps 


(5) 


The primary source of variation in touchdown point will be due to 
the effect of steady winds and thus in the following idealized analyses 
we consider only steady horizontal winds (independent of altitude.) 

a. Velocity Relations 


The basic velocity relations are shown in Fig. 12 according to the 
above assumptions. This diagram is shown for positive flight path 
angles, i.e., climb; however, the relations will also apply to descent 
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shown for 
X>0 ( climb) 



Vg 

Figure 12* Velocity Relations 

flight conditions. Since the flight path angle y will always be small 
(i.e., a few degrees) the ground speed may be readily approximated as 
simply the algebraic sum of the wind speed and aircraft airspeed 


Vg 

- 

<6a) 


= V„ + Va 

(6b) 

where is the velocity of the airmass with respect to 

is the velocity of the aircraft with respect to the moving 
sink rate is then given by Eq. 7 

the ground, 

a 

airmass. The 

• 

h = 

Vq tan Y 

(7a) 

« 

(Vw + Va)Y 

(7b) 


b. Horizontal Distance Traveled 


The horizontal distance traveled in the deceleration region is 
directly affected by windspeed and is 
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X 


+ 125.7 + 0.0634h)dt 




t 


= (Vw + 


125. 7)t + 0.0634 



(950 - 8.333t)dt 


= (Vv, + 185. 9)t - 0.2642 


( 8 ) 


Deceleration region trajectories computed from Eqs. 3 and 8 are 
plotted in Fig. 13 for the no wind condition and for 30 kt head and tail 
winds. The curvature seen in the trajectories is a consequence of the 
steady reduction in airspeed according to the schedule of Fig. 11 com- 
bined with a constant sink rate. The horizontal distance traveled in 
the deceleration region is found from Eq. 8 evaluated at t^ = 96 sec. 

Xo = 96 V„ + 15,412 ft (9) 

Xo is plotted in Fig. 14 as function of windspeed. The flight path 
angle at the end of the deceleration region, Yo is found from Eq. 7b and 
as 


To 


h ^ 8.333 

Vw + + 135.2 


( 10 ) 


and is also plotted as a function of windspeed in Fig. 14. 

c. Idealized System Response in Flare 

In the flare region the idealized sink rate response is (from 
Figs. 9 and 11): 
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“^BARo “ “ ^biaS 


= -h - 0.0505h - 0.75 = 0 

( 11 ) 

= h + 0.0505(h + 14.9) 

= h + i (h + hg) = 0 
where the flare time constant is: 


JL = 1 

0^0505 


19*80 sec 


and is set by the 80 kt airspeed, no wind flare condition. 

The ordinary first order differential equation in Eq. 11 forms an 
Initial value problem to define the flare trajectory when combined with 
two initial conditions. 


h(0) = ho = 150 ft (12) 

h(0+) = h(O-) = -8.333 fps (13) 

where now t is measured from the start of flare. It should be noted 
that while the sink rate in Eq. 13 is independent of wind the corre- 
sponding flight path angle will not be. These two initial conditions 
require the continuity of altitude and glideslope at the deceleration/ 
flare interface. Laplace transforming Eq. 11 


(sh(s) - h(0)) +4 h(s) = -4^ 

T T 


(14a) 


h(s) = 




(s + 1 /t) s(s + 1/t) 


( 14b) 
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Inverse Laplace transforming gives the altitude as function of time 


hoe-t/t - hg(l - e-t/x) 

(15a) 

(ho + hg)e-t/x - hg 

(15b) 

164.9 EXP(-t/19.8) - 14.9 

(15c) 


d. Time for Touchdown tj 

The elapsed time to touchdown measured from the start of flare, t^ 
is derived as follows. 


h(ti) = (ho + hg) e - hg = 0 


(16a) 


0 A Si - - - 

ho + hg 


(16b) 


tl 


-T)ln 



+ hg 


} 


-19.80 £n 


r 1 

450 + 14.9'' 


47.6 sec 


(17) 


e. Altitude Bias, hg 

The altitude bias, hg, in Eq. 11 is used to set the altitude below 
the runway for the zero sink rate asymptote (see Fig. 15) positive 
values of hg, insure that the sink rate at touchdown will be greater 
than 0, i.e., h(t^) < 0. 
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Figure 15* Definition of hg 


f* Sink Rate at Touchdown 


From differentiation of Eq. 15 the sink rate at touchdown is 


•/ . ^^o -ti/x 

h(ti) = - ( J e i 


_>o h. 

^ T- + hB^ 




= cons 


tan£ 


(18) 


From the preceding, it may be seen that the altitude variation, the time 
for touchdown, and the sink rate at touchdown are all independent of 
windspeed* 

g • Horizontal Travel in Flare 


The horizontal distance traveled during flare (measured from flare 
Initiation) is 
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X(t) 


(Vw + 135. 2)t 


(19) 


= Vet = (V„ + Va^)t = 


where Va^ is a constant 135.2 fps (80 kts) during flare. The distance 
traveled in flare to touchdown, is 


Xi = VGti 


= (V„ + 135.2K47.6) 


( 20 ) 


Thus, the distance traveled in flare is a function of windspeed and is 
plotted in Fig. 14. The total distance traveled in the deceleration and 
flare regions is obtained by summing Eqs. 9 and 20, i.e. 


Xt - Xo + Xi 

= 143.6 + 21,848 


( 21 ) 


The total distance is also plotted in Fig. 14. 
h. Summary 


While this mechanization produces quite acceptable sink rates at 
touchdown, the landing dispersions are highly dependent upon initial 
approach alignment and the direction and magnitude of atmospheric dis- 
turbances. Furthermore, the ground pilots view of the landing through 
the TV monitor will vary from landing to landing because of these same 
alignment and atmospheric factors. This makes it difficult for the 
pilot to monitor for proper autoflare and generally contributes to the 


TR-1 180-1 


29 



requirement for several practice approach and landings with the PA-30 
aircraft just prior to each HIMAT flight. 

2. Lateral/Dlrectlonal Hechanlzatloa 

PA-30 autoland roll and yaw axis mechanization block diagrams are 
presented in Figs. 16 and 17, respectively. 

Unfortunately, there is no onboard lateral flight path sensor for 
the roll or yaw axis during autoflare as is the case for the longitudi- 
nal axis. Therefore the ground pilot must continue to perform this 
function via the T.V. display and uplink commands to the vehicle. The 
uplink inputs are discrete left or right yaw rate commands. The HIMAT 
backup control system is based on body angular rate sensors and has no 
direct measure of bank angle. However, since turn rate is directly pro- 
portional to bank angle in a coordinated turn, the yaw rate may be sub- 
stituted for roll attitude. The roll axis mechanization is thus a turn 
rate command, zero turn rate hold system. 

The ground pilot inputs discrete turn rate commands until the 
desired directional path is achieved and then, with the command removed, 
the roll axis controller attempts to maintain zero yaw rate from that 
desired directional path. The integral of yaw rate error is fed back to 
provide a tight control structure at low frequency, l.e., trim. 

As noted previously, the yaw axis (Fig. 17) consists of a conven- 
tional washed out yaw rate damping feedback and a turn coordination cir- 
cuit which combines washed out yaw rate with lagged lateral accelera- 
tion. That this provides turn coordination may be seen as follows: 

Ki _ K2S Ki ^ K2 

(s + 1) (s + 1) “ Js + 1) W kJ" ^ J 


Now if the lateral acceleration sensor is located along the vehicle 
x-axis and at the center of rotation for forces applied at the rudder, 
then it senses only vehicle translation and not rotation due to rudder 
deflection. This is the ideal sensor location for the signal to be used 
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Figure 16. Roll Autoland (From Reference 6) 



Figure 17. Yaw Autoland (From Reference 6) 
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for turn coordination and to be uncontarainated by rudder induced yaw 
acceleration due to the damper function. In simple terms 


aycg = Yyv + Y5j.6r 

^ys ° ^ycg 

= Yvv + Yg^fir + Xar 


and the accelerometer is located at a distance from the c.e. for 
which the X^r component will tend to offset the component. Then 

one is left with YyV (or Y^3) the major part of the acceleration 
sensed. If the sensor cannot be placed at this ideal location due to 
physical constraints, then an equivalent signal can be synthesized by 
summing the sensor signals as noted above and with K 2 /Kj^ = X^. The 
resulting signal is then integrated to to maintain zero v (or g) trim. 

This overall lateral/directional mechanization would probably be 
adequate under absolutely calm wind and full, three dimensional vision 
conditions. However, in the presence of turbulence, crosswind, or wind 
shear, the bang-bang nature of the command makes precision path control 
impossible. The task is further complicated by the two dimensional dis- 
play characteristics noted previously in Subsection A. 3. 

The Ref. 6 early PA-30 flight checkout verifies this assessment by 
specifically commenting that 

"it was difficult to control roll attitude, combat 
drift, etc." 

and 

"the major piloting task of the autoland mode is the 
lateral heading control task." 

An example of the lateral control imprecision and pilots' workload is 
given in Fig. 18 which has been taken from the Ref. 6 flight records. 
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SECTION III 


IMPROVED AUTOLAND CAPABILITY 


One goal of this study is to synthesize an "improved" terminal con- 
trol autoland capability for RPRV^s (and the PA-30). The improvements 
sought involve both consistency of terminal conditions and easing of 
ground pilot workload. Further, in the interest of timely implementa- 
tion potential and overall cost, it is desirable to make as much use of 
existing DFRF equipment and facilities as possible. 

The essence of any good performing automatic landing and flare sys- 
tem lies in making the aircraft tightly track an approach beam and then 
wind proofing the flare to a smooth landing at a consistent touchdown 
point and touchdown sink rate. Thus, two separate automatic modes are 
involved in the "improved" autoland system: automatic beam tracking to 
produce consistent flare initiation conditions and an automatic flare 
mode which minimizes the influence of head, tail, and crosswind condi- 
tions encountered during the approach and flare. An additional consid- 
eration is minimization of any transients due to switching from the beam 
tracking to the autoflare mode. The end result is a vehicle/automatic 
controller system which greatly reduces ground pilot visual and mental 
workload by providing essentially invariant approach and landing condi- 
tions to monitor. 

For the automatic beam tracking mode it is proposed that the present 
ground pilot station synthetic ILS flight director concept be extended 
to use the flight director needle command deflection signals as the 
basic deviation (or error) input to the flight control mechanization 
hosted in the ground computer. Otherwise the system operates on down- 
linked signals from onboard sensors and up-linked surface position com- 
mands as is now accomplished. Figure 19 pictures the concept in simpli- 
fied block form consistent with Fig. 1. 

For the Improved autoflare mode the "variable t" algorithm and mech- 
anization developed by Boeing (Ref. 7) for the B-737 aircraft is 
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HIMAT 

or 



Modified HIMAT BCS : add 0 ILS 

0 Improved autoflare 

Modified HUD : add flight director and status 

information 

Figure 19. Improved Autoland System 
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adopted. In this system the time constant of the flare law (x) is 
varied with ground speed to produce essentially invariant touchdown dis- 
tance and sink rate from a given flare initiation height despite the 
presence of head or tail winds. 

One of the main sources of touchdown dispersions for conventional 

* 

flare laws using pre-fixed h and h gains is the sensitivity to approach 
ground speed as affected by steady winds. This sensitivity is due to 
the time exponential nature of altitude and sink rate decay during 
flare, resulting in a fixed flare duration from a given altitude. The 
flare distance from a given altitude will therefore be proportional to 
ground speed, which in turn depends upon the wind. A high flare initia- 
tion altitude compounts the problems. 

Variation of h or h gains with ground speed produces different time 
exponential (variable x) decay of altitude and sink rate such that lon- 
gitudinal flare distance remains fixed at a cost of slight variation in 
sink rate at touchdown. 

Incorporation of ground speed correction generally is considered to 
require an onboard inertial navigation sensing package which would be 
too costly and bulky for RPRV (or PA-30) application. Fortunately, 
Boeing (Ref. 7) has already addressed this problem for their B-737-200 
and developed an autoflare law which utilizes a "pseudo” measure of 
ground speed. Additional features of the Boeing longitudinal mechaniza- 
tion are the provision of pre-flare glideslope tracking, transient free 

transition from glideslope tracking to flare initiation and incorpora- 

. 

tion of runway referencee, radar altitude rate (hp^) for the flare. 
Thus, this mechanization incorporates all the features desired for the 
RPRV improved longitudinal approach and autoland modes. For the lateral 
axis a rather simple beam tracking, wings level, flare mechanization 
such as employed in the Lockheed L-1011 and Compass Cope is proposed. 

The following subsection describes the longitudinal mechanization in 
more detail. Idealized system performance is then presented in the next 
subsection for comparison with the current autoland responses discussed 
in Section II. The third subsection describes the proposed lateral sys- 
tem mechanization. 
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A. LONGITUDINAL BEAM TRACKING AND ADTOFLARE 


1* Mechaoizatioa 

A block diagram of the proposed mechanization is presented in 
Fig, 20. Except for the glideslope error (h^gg) input, all sensed 
parameters are the same as in the present autoflare mechanization 
(Fig. 9). Switching and signal fading between modes is on the basis of 
a pre-selected, runway referenced, flare initiation radar altitude, hj^^, 
as follows: 

a. hR > 

The beam tracking mode is a conventional proportional plus integral 
system. Note that when hj^ > K(^SE = ^RL constant (limit 

value) and the integral gain is integral path nulls the 

beam deviation error while the lagged normal acceleration and pitch rate 
feedback paths provide, respectively, heave error and short period atti- 
tude-damping. The latter two feedbacks are of the same form as in the 
present autoland deceleration and flare mechanization. 

. 

The additional derived sink rate feedback path, h, provides the 
"pseudo” measure of ground speed. This may be noted by observing that 
in the steady state (above 

hRL = constant 

hRL (tps'%'1^ = ° 

• • ( 1 ) . 

^ = ^bias + j (Tps + T) + (Tps -h 1) ) ^CF 

, (Tps + 1) , 

= hbias + 1) ^CF 
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+ above 



Figure 20. Integrated Glide Slope/Variable Tau Flare Control Law 








hbias + hcF 


Thus, h is proportional to complementary filtered barometric altitude 
rate as the vehicle descends the glideslope plus a small bias to assure 
a firm touchdovm rate* Since the glideslope is fixed with respect to 

the ground, the rate of descent along the glideslope is proportional to 

• • 

vehicle ground speed, i.e. h = Vq tan y* The h feedback provides a beam 
deviation bias which initially causes the vehicle to fly above the 
glideslope until the h signal is cancelled by the integral of beam error 
path. Therefore the output of the integrator becomes proportional to 
vehicle ground speed in the steady state descent along the glideslope. 

b) hj^ < h^ — When the flare initiation altitude is reached several 
things happen simultaneously to rapidly and smoothly transition the sys- 
tem into the autoflare mode. First, in the beam tracking path, the 
gain K^sE rapidly faded to zero so that the integrator output, Kj^/s, 
holds its output (which is proportional to ground speed, V^, at the time 
of changeover). Secondly, the hj^^ signal comes off its limit value and 

starts to decrease. Thirdly, the switch opens in the first order lag 

. 

path of the complementary filtered altitude rate, h^>p, and the signal in 
this path starts to decay at a time constant T^, At the same time the 

decreasing radar altitude, h^j^, passing through the washout circuit pro- 

. 

duces a signal proportional to hj^. The two rate signals complement one 
another to produce up elevator to initiate the flare, i.e.. 


t- 

Tf 2Tf 3Tf 
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Thus, the signals produce a smooth transition to an autoflare system 
that is completely runway referenced within three time constants, Tp, 
after flare initiation. 

The variable x autoflare block diagram may then be considered in the 
simplified form shown in Fig. 21. The elevator command for flare, 

Is 


6 




TpS 


+ -(ximr TTf's + 1) + hfilAS 


=■ + Hr + TphcF + hsiAS 


- ^hR^G ^ ^BIAS + TphcF 


= ^hj^VcCT s + 1) hR + hsiAS + TphcF 


Vg 



Figure 21. Simplified Autoflare Block Diagram 
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where x = is the ground speed ^variable tau" and the contribu- 

tion through the complementary filter is primarily a smoothing term. 

The end result is a full runway reference, ground speed adjusted 
flare law. Example flare trajectories from the Ref. 7 B-737 performance 
assessment with wind variation from -15 to +25 kts are shown in Fig. 22. 
While the insensitivity of the touchdown point to wind is quite appar- 
ent, it should be noted that the slight variation in flight paths will 
result in similar variation in sink rates at touchdown. These aspects 
will be assessed in greater detail in the next subsection. 

As a final note, in addition to achieving the desired performance 
goals the mechanization also is compatible with the existing HIMAT BCS 
(and PA-30) autoflare mechanization. Figure 23 shows the present auto- 
flare (mechanized for the PA-30) as solid lines and the paths and blocks 
necessary to incorporate the integrated glideslope/variable tau flare 
laws in dashed lines. The improved system requires only a minimum of 
added circuitry plus incorporation of the glldeslope error input. The 
circuitry could be built into the onboard backup control system with the 
GSE signal transmitted via up-link. Since all other sensors are already 
onboard, loss of the up/down communication link could automatically 
result in the system reverting to the present HIMAT/PA-30 autoflare con- 
figuration as a part of an autononous landing configuration. 

2m Idealized Performance 

The present system is subject to significant variation in touchdown 
point in steady winds as shown in Figs. 13 and 14. In this section it 
will be shown that this problem may be essentially eliminated by use of 
the proposed system which employs automatic tracking of the synthetic 
ILS followed by a "variable x" flare law. 

The idealized trajectory for the proposed system is shown in 
Fig. 24. It will be demonstrated that for the proposed system: 

• The trajectory is completely defined by specify- 
ing Yo> airspeed (V^ ), zero wind flare 

time constant (Xq) and touchdown sink rate 
(h(ti)). 
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GEAR HEIGHT VS DISTANCE DURING FLARE, 
WINDS : 
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Figure 23. Compatibility of Present and Proposed Autoflare Mechanizations 










Figure 24. Trajectory for Proposed System 

9 Once defined, the trajectory shape is invariant 

with wind, in particular, h^ is independant of 
wind but does vary with choice of Xoj V;, _ , Yn> 
and h(ti). ° 

• Unlike the existing system, touchdown time and 
touchdown sink rate do vary with wind velocity, 
but not significantly* 

• The proposed and existing system always have 
(slightly) different flight paths before flare if 
deceleration occures during the approach phase* 


Synthetic ILS Beam Following 


For purposes of comparison with the preceeding analysis of the 
existing system, it will be assumed in the following that the ILS track- 
ing region coincides with the deceleration region and that transition to 
flare is made at h^ = 150 ft as in the existing schedule* These assump- 
tions are not necessary and perhaps not even desirable for the actual 
implementation of this system but will be used here only for purposes of 
comparison* It is shown in Appendix D that the speed and altitude 
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schedules can generally be considered independent and that the beam 
tracking system should perform adequately either at constant speed or 
with a scheduled speed variation in the region between 110 kts and 
80 kts as used in the present PA-30 schedule, Fig. 11. The total dis- 
tance traveled in the beam tracking deceleration region, X^, is: 


Xo 



( 22 ) 


where h^ = h^ (tq, Yq, h(t]^)) will be determined in the next sec- 

tion. 


b. Variable x Flare Law Ideal Altitude Response 

As noted previously, the 'Variable x” flare law concept is used to 
fix the flare trajectory in space essentially independent of wind by 
making the flare time constant x a function of windspeed. The time con- 
stant X is therefore redefined as 



(23) 


where Xq, Vq^ are the no wind parameters. Following Ref. 8 Eq. 23 is 
substituted into Eq. 11 multiplied by l/K^^ giving 


^oVg, 


° + (h + he) = 


Vg dt 


(24) 


Since 


dh _ ^ ^ 

dt dx dt dx ^ 


(25) 


Eq. 24 becomes 


*^b) = 0 (26) 

Equation 26 is an ordinary differential equation in spatial variables 
for the flare trajectory and is independent of steady wind. By analogy 
to Eq. 15b 
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h(x) 


(27) 


= (ho + hg) 


-x/toVg 

e 




where x is measured from the flare initiation point, 
c. Initial and Final Trajectory Points 

Since Eq. 27 is independent of wind, the flare initiation height, h^ 
and touchdown distance, measured from flare initiation are indepen- 

dent from wind also. From Eq. 24, evaluated at the start of flare, 

ToVGo . 

ho — ho - hg 

^oVGo 
’ G 




(28) 


(recall that = Vq^) 
Thus, Eq. 22 becomes 



From Eq. 27 evaluated at touchdown (h = 0) 

■ (to va“ Yo) 


(29) 


(30) 
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80 kts 


Using the PA-30 values of Tq = 19.8 spec, hg = 14.9 ft, and Vag = 
as in Section II; X^, and X^ are plotted as functions of the 

glideslope at flare initiation Yq in Fig* 25. The ideal flare trajec- 
tory as defined by Eq. 27 is identical for all wind speeds, however, the 
actual system with additional feedbacks and finite bandwidth for the 
various feedback loops will show some slight variation in trajectory 
with windspeed, e.g.. Fig. 22; however, this variation is very small 
compared to that for the existing BCS system and thus the preceding 
analysis is quite adequate for assessing the proposed system. 

d. Comparison Trajectories for Proposed and Existing Systems 

It is difficult to directly compare the proposed system with the 
existing BCS autoflare as implemented on the PA-30 because even the 
idealized trajectories have slightly different shapes. Specifically the 
constant h region in the existing system produces a curved trajectory if 
deceleration occurs (see Fig. 13). Probably the most meaningful way of 
comparing the two systems is to specify the same sink rate at touchdown 
(no wind) for both systems. As will be shown shortly, this will occur 
if h 3 and x = Tq are the same for the two systems. If the flare air- 
speed, and the flight path angle at flare initiation are also the 

same, the shape of the flare trajectory of the proposed system in any 
wind will be the same as for the existing system with no wind. This may 
be seen in Fig. 25 by noting that the h^ and values for the proposed 
system correspond to the existing system values with no wind at the 
appropriate Yq -3.5 deg. 

Because of the slight difference in shape of the trajectories before 
flare in the two systems, they differ in and Xp for Yq = “3.5 deg. 
It may be seen that the total travel for the proposed system will be 
several thousand feet less. 

e. Sink Rate at Touchdown for the Variable x System 

As noted above, spatial quantities such as touchdown distance and 
flare initiation height are ideally independent of wind speed for the 
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variable x system. However, time dependent quantities such as sink rate 
at touchdown and touchdown time do vary with wind speed, a situation 
which differs from the present BCS autoflare law. The sink rate at 
touchdown may be derived from Eq. 24 by setting h = 0. 


h(ti) 


^bVg 

^oVGo 



(31) 


As may be seen from Eq. 31, the sink rate at touchdown varies directly 
with windspeed. Since the slope of flight path at touchdown is indepen- 
dent of wind speed, the difference in sink rate is due to the direct 
variation of ground speed with wind speed for a constant airspeed. The 
effect is not significant as shown in Fig. 26. 

HEADWIND Horizontal Wind Speed tailwind 



Figure 26. Variation of Touchdown Sink Rate With Wind 
for the Proposed System 
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f. Touchdown Time for the Variable t System 


The touchdown time for the variable x system Is computed as shown 
below. 

Xi = X(ti) = (Vw + Vao) ti 

From Eq. 29 

ti = Xi/(V„ + Vao) 

-Tq Vap \ / -hB 

+ Vsq J \ ^o ^®o Yo 

<0 <0 

As shown In Eq. 32, touchdown time varies Inversely with wlndspeed such 
that the airplane lands sooner with the tall wind, as the consequence of 
the higher groundspeed and a constant Xj^. This Is In distinction to the 
present BCS autoflare system where the elapsed time In flare Is 47.6 sec 
Independent of wind. 

Thus, In summary, the primary effect of the beam follower and the 
variable x flare law Is to trade the temporal Invariance of the existing 
system for spatial Invariance. The value of this tradeoff is that the 
aircraft will always land at essentially the same point on the runway. 
As a consequence, the same landmarks may always be used by the ground 
pilot to monitor correct function of the approach and autoflare systems. 

B. LATERAL BEAM TRACKING AND FLARE CONTROL 

Lateral control is considerably simpler than longitudinal In that 
synthetic beam (runway centerline) lateral deviation information can be 
derived from the radar tracking signal and ground computer and transmit- 
ted to the flight control system throughout flare as well as approach. 
However, the deviation accuracy and signal noise may degrade somewhat as 
the vehicle nears the runway. Thus, the vehicle can be automatically 
controlled to the runway centerline until very close to touchdown and 
therefore minimize lateral drift due to crosswinds. 



> 0 


(32) 
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Past PPRV's (HIMAT, 3/8 scale F-15, etc.) have generally employed 
skid type landing gear which is designed to be sufficiently sturdy to 
survive crabbed landings and thus avoid complicating the manual control 
task just prior to touchdown. There appears to be no reason to alter 
this landing technique and unnecessarily complicate the flight control 
mechanization. 

Thus, a simple, proven lateral mechanization such as used in the 
L-1011 (Ref. 9) and Compass Cope RPV (Ref. 10) is proposed. This mecha- 
nization, Fig. 27, employs beam error, lateral acceleration, and roll 
attitude inputs to aileron control. There are no crossfeeds to rudder. 
The beam error is integrated to minimize effects due to crosswind. Rate 
of beam error is generated by passing the signal through a washout cir- 
cuit. The rate signal is summed with lagged lateral acceleration (from 
a body mounted accelerometer) and roll attitude to produce a complemen- 
tary filtered damping term. This complementary filter may not be 
required for the RPRV facility however. It is noted in Ref. 9 that if 
the ILS beam is fairly clean the complementary filtering is not neces- 
sary. Since the "synthetic ILS" beam is computer generated and trans- 
mitted to the flight control system at a known sample rate, "beam noise" 
should be quite narrow band and readily handled. 

After flare initiation the beam error gain is reduced to zero as a 
function of radar altitude (runway referenced) so the roll attitude con- 
trol loop becomes dominant as the vehicle approaches touchdown. The 
roll command becomes zero and assures a wings level touchdown so that 
there will be no danger of a wing strike. As with the longitudinal 
axis, the vehicle lateral response during flare will be consistent in 
each landing and therefore reduce the pilot monitoring workload. Of 
course, lateral drift due to cross winds will build during the wing 
leveling phase and will be dependant upon specific conditions of each 
landing. 

It appears the only way to assure a wing level touchdown is to feed- 
back bank angle from an airborne vertical gyro. This should be no par- 
ticular problem since the HIMAT does carry a V.G. as a part of the emer- 
gency recovery system (but its output is not used for active control) 
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Figure 27. Proposed Lateral Autoland 
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and the PA-30 apparently employs a V.G. roll output for the manual roll 
command controller (Fig. 7 and Table 1). Presumably any future RPRV^s 
'^lll have similar needs for attitude signals. It is only necessary that 
the sensor output be employed in a direct manner in the autoland beam 
tracking and flare modes. 

Integration of the improved lateral autoland control, a manual roll 
rate command mode, and the present HIMAT backup control system (BCS) is 
shown in Fig. 28. The manual roll rate command mode is employed to fly 
the vehicle to localizer beam acquisition. The autoland mode is then 
engaged to track the beam through flare. In the event of uplink or 
downlink failure, the system will revert to the present backup control 
mode with turn commands, r^, from either the current chase plane control 
mode or to the autonomous navigation and orbit mode discussed later. 

The rudder axis control system can remain the same as in the present 
autoland system since its function is to provide yaw rate damping and 
turn coordination . 
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Figure 28. Possible Integration of Automatic and 
Manual Lateral Control Modes 
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SECTION IV 


PROPOSED HEAD-UP APPROACH AND LANDING DISPLAY 

The current (Ref* 4) head— down ADI and head— up T*V. displays and 
their shortcomings were summarized In Section II- In this section a 
relatively simple scheme Is presented to enhance the head-up T.V. dis- 
play to alleviate current problems. The specific goals of the Improved 
display are to 

• integrate appropriate cockpit panel data onto the 
T.V. display so that the pilot can fly the com- 
plete approach and landing using the head-up dis- 
play as the primary visual information. 

• Incorporate status and flight director informa- 
tion to minimize pilot workload in monitoring 
automatic beam following and landing flare system 
performance and in accomplishing the task 
manually. 

Means of enhancing T.V. display for the accomplishment of these 
goals have also been addressed in Refs. 11 and 12. For example. Ref. 11 
describes advantages and limitations of several different glide slope 
type Information aids superimposed on the T.V. display from body-fixed 
versus Inertially stabilized camera mountings. The display aids 
included fixed reticle, glide slope reference lines augmented by rate 
and acceleeration Information. The latter is in effect a beam director 
display. A piloted simulation involving three different sizes and types 
of aircraft showed the beam director to result in the least displace- 
ment, attitude, and control force variances with either the body-fixed 
or inertially stabilized camera systems. The director algorithms were 
of the form 


^GS = hg + a^hg + biag 

£LOC = 7e + ^lye + 
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where and are the director glide slope and localizer devia- 
tions, respectively. As will be shown later, this algorithm does not 
include all of the feedbacks considered necessary in Ref. 13 for good 
pilot-centered display but it is compatible with the automatic beam fol- 
lowing mechanizations proposed in the previous section and tends to 
indicate a simple display of the elevator command would probably ease 
the pilot workload during the approach - 

Reference 11 did not address the landing flare maneuver. However, 
this aspect has been quite adequately covered by the Ref. 12 simulation 
investigation in which several parameters were superimposed on the T.V. 
landing view. This investigation even used the PA-30 aircraft as the 
example vehicle. It was determined that the most consistent flare and 
touchdown was achieved with a longitudinal dynamic display of radar 
altitude and altitude rate. The display parameters were ratioed in much 
the same manner as the two major feedback paths of the automatic flare 
law previously presented in Fig 21. 

The Ref. 12 simulation also investigated the usefullness of glide 
slope and localizer information during the approach phase and came to 
much the same conclusions found in Ref. 11. However, Ref. 12 incorpora- 
ted various flight condition and attitude cues to assess the added bene- 
fit (or confusion) aspects. It was found that inclusion of airspeed, 
altitude, and altitude rate digital readouts were of definite benefit. 
Recommendations were made for the location of these parameters on the 
T.V. monitor. A horizon line and a roll attitude indicator were sepa- 
rately assessed. Attitude information is generally considered to be 
highly beneficial by providing lead on path control and thereby reducing 
pilot workload while improving performance. It was found in Ref. 12, 
however, the roll attitude indication was ignored by the pilots. This 
may have been due to the particular display symbology employed and the 
visual clutter it introduced. It was determined that inclusion of the 
horizon line was beneficial. 

Based upon the results of Refs. 11 and 12, the approach director 
display format of Fig. 29 has been selected. In this example the air- 
craft is below the synthetic glide slope beam (^ 0 ^ 3 ) and to the left of 
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the runway centerline (FD^^q^,). The horizon bar indicates the aircraft 
Is pitched down and banked left. To prevent the horizon bar from ever 
leaving the T.V. screen or interfering with the other display param- 
eters, the pitch gain should be very low. The horizon bar thus becomes 
primarily a display of bank angle. Airspeed, altitude, and altitude 
rate are located as recommended in Ref. 12. Figure 30 is an example of 
the approach display augmentation superimposed on the different T.V. 
views of Fig. 5. These assume the aircraft to be in the same location 
with respect to the runway in each case but with different bank and 
heading angles. By flying the vehicle reference cross to the flight 
director bars the vehicle will follow the synthetic ILS beam to the run- 
way. 

The pilot should be alerted to approach to altitude for flare initi- 
ation by a flashing symbol. This might be the vehicle reference cross 
or the altitude numeric. A few feet above flare initiation altitude the 
display should change to the format indicated in Fig. 31. Here the run- 
way localizer director bar and airspeed numeric continue the same func- 
tion as in the approach mode. The horizon bar is removed and replaced 
by a runway reference and the two glide slope bars move independently of 
one another with the right hand bar indicating radar altimeter altitude 
from the runway and the left hand bar Indicating rate of descent. Bar 
deflection gains are s.elected so that the desired exponential flare 
results in the bars being aligned (horizontally) and simultaneously 
approaching and reaching the runway bar at vehicle touchdown. The 
example of Fig. 31 indicates the sink rate to be excessive for the 
immediate altitude. Log scales are used to Increase display sensitivity 
as the vehicle approaches ground level. 

The following subsections address in more detail Implementation of 
flight director beam following and landing flare algorithms and an opti- 
cal means of superimposing the video and computer generated information. 

A. BEAM FOLLOWING MONITOR/DIRECTOR MECHANIZATION 

In the automatic beam following mode the ground pilot must monitor 
for proper error nulling and corrective action for any large deviations 
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due to external disturbances. This function can be served by display- 
ing, for example, an indicator deflection proportional to beam error 
alone. However, in the manual mode the display must take into account 
various pilot-centered director requirements including 

• ability of pilot to distinguish between responses 
due to external disturbances and his corrective 
inputs . 

• elimination of any integral of beam error feed- 
back so the display does not ramp off to its 
limits during any periods of unattended opera- 
tion. 

• provision of effective controlled element dynam- 
ics which approximate K/s in the bandwidth region 
over which the pilot is attempting control. 

The first two of these requirements define a display algorithm dif~ 
ferent from that employed by the automatic beam follower mode previously 
presented. To avoid any unnecessary and confusing transients when 
switching from auto to manual, or vice versa, it is desireable that one 
display concept be used for both modes. Thus, the manual mode should 
prevail. 

The basic considerations for selection of longitudinal flight direc- 
tor system feedbacks are summarized in Table 3 from Ref. 13. It should 
be noted that the key beam feedbacks are beam deviation and deviation 
rate (in our case hQ5g and hQgg). Beam error Integration is not desired 
as noted previously. Washed out pitch attitude is desired to provide 
short term feedback of vehicle response to pilot inputs and therefore 
prevent any tendency to over control. These three feedbacks provide 

precise beam following capability. Including pitch rate allows indepen- 
dent control of the short period damping ratio. 

A resulting beam-derivative plus pitch attitude system is shown in 
generic form in Fig. 32 where and are the respective feedback 
functions . 

A central consideration of the pilot-centered requirements is that 
the feedback gains and equalizations be selected so that the net 
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TABLE 3. CONSIDERATIONS FOR SELECTION OF LONGITUDINAL SYSTEM FEEDBACKS 
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Figure 32. Simplified Director/Vehicle Block Diagram 


transfer characteristics from elevator, 6^, to director instrument dis- 
placement, FD, look approximately like an integration, K/s. As shown in 
Ref. 13, the analytical procedure for assessing and establishing this 
result is to add the component vehicle transfer function with their 
associated equalization. The effective flight director transfer func- 
tion is then. 


FD 


^FD 0 FD h 

^0 + Gji Ng 


e 


A 


where and are the vehicle pitch and altitude response numera- 
tors, respectively- In keeping with the feedbacks of Table 3, the feed- 
back transfer characteristics are: 


,FD Kqs _ _ K§s[s + (l/T„o + Kq/K§)] 

® Is + U/T„oJJ “ (s + (l/T„o)J 


Gh = Kh + KhS = (s + (Kh/Kh)] 
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Figure 33 is a Bode-Siggie amplitude sketch of the director transfer 

FD 0 FD h 

function numerator, G 0 particularized for the PA-30 

dynamic and the Ref. 13 recommended equalization values (see Table 4). 

The numerator roots are determined by the intersection of the gain 

closure line and the various a and jw amplitude loci in Fig. 33. 

Selecting a value of +0.035 rad/ft for the gain K^/Kq produces effective 

flight director zeros identified by the symbols (§) and results in the 

FD 

flight director numerator shown. 

FD 

A sketch of the open-loop director/vehicle, Ng^/A, Bode asymptote 
and phase is given in Fig. 34. The selected director equalization and 
gains are seen to result in the desired K/s-llke amplitude ratio over a 
large frequency region where the pilot might be expected to close the 
loop (potential crossover region). The phase sketch for this linear 
analysis indicates there is no instability problem regardless of how 
high a crossover the pilot might try to achieve. Actually, however, the 
effective time lags (t) due to the up/down link communication can be 
expected to produce a phase curve such as the dashed line in Fig. 34. 
This still allows a bandwidth of 2 rad/sec without any pilot generated 
lead and somewhat greater than 2 rad/sec if pilot lead is employed. 

It should be emphasized that the parameters of this example longi- 
tudinal flight director law are particularized to the PA-30 dynamic 
characterized. For such a director to work properly with any other 
vehicle, the feedback gains, time constants, etc., must be rescheduled 
to match that vehicle. This is not difficult in that it merely requires 
repeating the foregoing analysis to select the gains and equalization. 
These are then set into the display algorithm in the ground computer for 
each vehicle. 

The lateral flight director algorithm is of the same form as the 
longitudinal except that the feedbacks are localizer ye and yg, washed 
out roll attitude and roll rate. These feedbacks are all consistent 
with the autoland mechanization of Fig. 27, however, the appropriate 
gains, time constants, etc., must be tailored to each vehicle using the 
same analytic approach discussed above. 
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(rod/sec) 
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Figure 33. Bode-Siggie Amplitude Sketch of the Director Transfer Function Numerator 








TABLE 4. SELECTED EQUALIZATION VALUES 


EQUALIZATION 

EXPRESSION 

DESIRED LOCATION 
AND VEHICLE VALUE 
(rad/ sec) 

SELECTED 

EQUALIZATION 

REMARKS 

Pitch 

Washout 



1 

^WO 

^ ^sp “ 3.51 
> ^ = 1.93 

T02 

2.5 

Washout less than 
^sp provide 

attitude stability 
but greater than 
1/^02 for wind- 
proofing and to 
maintain altitude 
bandwidth. 

Pitch 

Attitude j 

i 

I 

i 

! 

i 

1 ^ Kq 

Two K0 

^ ^sp “ ^ 

4.0 

Pitch attitude lead 
set to improve 
short period damp- 
ing and extend the 
region of K/s by 
having the result- 
ing (00 zeros cancel 
the o)gp poles 

Altitude i 

! 

Kh 

! 

1 

= o)p = 0.16 

0.2 

Greater than (Op to 
avoid a "busy" dis- 
play and the low 
frequency closed- 
loop d/d^ amplitude 
droop, yet maintain 
mid-frequency phase 
margin 
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PA-30 Aircraft Dynamics 





B. LANDING FLARE MONITOR/DIRECTOR MECHANIZATION 


The landing flare is basically a terminal control problem to main- 
tain the proper relationship between h and h to touchdown. Therefore, 
for monitoring the longitudinal autoflare system (Fig* 21) the two com- 
ponents of the elevator command, i.e., K^j^VQh^ and h, are displayed as 
in Fig. 31. This provides lead Information to the pilot in that he can 
independently discern sink rate and altitude above the runway. If the 
two signals simultaneously approach zero the autoflare is functioning 
properly. If sink rate should become excessive for any reason, the 
pilot can immediately detect it, take over manual control, and correct 
the situation. 

For manual flare control, there again should be no change in the 
display format from that employed in the automatic mode. Fortunately, 
there is no problem in that the same h^ and h parameters are consistent 
with the manual display recommended from the Ref. 12 simulation. Thus 
upon reaching the flare initiation altitude the longitudinal monitor/ 
director display signals are obtained from the Fig. 20 autoland block 
diagram. This is shown in simplified form in Fig. 35. 


Vg 



Figure 35. Landing Flare Monitor/Director Mechanization 
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The lateral axis control task during flare is to maintain the air- 
craft flight path along the runway centerline. This is accomplished by 
continuing tracking to the synthetic ILS beam all the way to touchdown 
(or to the lowest altitude for which the "beam" can be adequately gener- 
ated). Thus, the lateral beam tracking function remains the same for 
flare as for approach and monitor/director mechanization is unchanged. 

C. DISPLAY IMPLEMENTATION 

There are several possible means of combining the monitor/director 
displays with the present T.V. monitor. One is to electronically super- 
impose the two signals directly into the monitor screen. A second, less 
expensive approach is to optically superimpose the T.V. scene and a CRT 
generated display via a combining glass. This latter has been employed 
very successfully in several simulators fabricated at STI. This 
approach does require a reduced light level in the room containing the 
simulator. It is understood that the ground pilot station is presently 
operated in semi-dark conditions for better T.V. viewing and to minimize 
distractions in the pilots^ peripheral view. Thus, the combining glass 
approach is recommended. 

A possible layout is sketched in Fig. 36. This requires only the 
addition of a CRT and combining glass to the present ground station. 



Figure 36. Proposed Means of Superimposing Flight 
Director Symbology on the T.V. View 
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The combining glass is set at a 45 deg angle in front of the T.V. moni- 
tor. The display symbology of Figs. 29 and 31 is generated on the CRT 
and projected onto the combining glass. The pilot then sees the two 
views superimposed. This approach offers a distinct advantage in that 
the CRT display brightness can be independently controlled to optimize 
contrast between the real world T.V. scene and the director symbology. 
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SECTION V 


AUTONOMOUS ORBITING ALGORITHMS FOR RPV^S 


A. INTRODUCTION 


This working paper documents the derivation of simple, orbit-about- 
a-f Ixed-point algorithms for autonomous navigation of a remotely piloted 
vehicle (RPV). This can form a preliminary basis for one of the on- 
board navigation software modules. 

The basic problem is one of providing a "safe orbit" capability in 
the event of loss of ground-to-air communications links. The autonomous 
mission phases or "modes" consist of: detecting uplink loss and deter- 

mining current vehicle and wind states; proceeding to the safe orbit 
area; initiating and maintaining an orbit about a fixed ground point; 
and after a predetermined number of orbits, exiting the orbit mode and 
proceeding toward a possible autonomous landing. Of course, any of the 
autonomous modes could be interrupted in the event of uplink recovery, 
and should allow for return to manual control. This working paper dis- 
cusses the orbit maintenance mode only. 

The constraints assumed are that: 

• The vehicle^ is subjected to winds aloft, which 
without correction would cause drift out of the 
safe orbit area, during the 5 to 15 minute orbit 
mode period. 

• The only information available to the vehicle is 
its position and an estimate of the wind vector 
at uplink loss, and continuous signals for air- 
speed and Euler heading angle (from onboard gyros 
or integrated rate gyros). 

• The module output is roll angle command (of 
course this is kinematically similar to a yaw or 
heading rate command). For simplicity — and in 
view of the "emergency" nature of the mission — 
speed control via autothrottle is not considered, 
and airspeed and altitude, for the moment, are 
assumed constant. 
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For this preliminary analysis, sensor noise and drift, and changes in 
the wind vector are assumed to be zero. The vehicle model is kinematic 
(with no sideslip or dynamics). Further analyses should look at the 

effects of these variables on orbiting accuracy. 

The remainder of this paper describes two algorithms. The first is 
for a circular path; the second, improved one is for dead reckoned wing 
pointing, which results in superelliptical orbits. Both are given in 
analytical terms, unembellished. Some trajectory plots are given, 
resulting from digital simulation. The appendices give time histories. 

B. CIRCULAR ORBIT ALGORITHM 

1* Problem: Fly Circular Path Over Ground » Under 

Steady Wind 

Assume : 

• Constant airspeed magnitude, 

• Zero sideslip, no vehicle dynamics 

• Heading angle, available 

• Orbit begins at t=0'*’ 

• Initial flight path is tangent to desired circle 

• Constant (inertial) wind vector (V^) 

Now, referring to Fig. 37, the ground speed is. 
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X = 


cos 


-1 (U + V cos \|;) 

Aj2 + 2UV cos t + V2 


the (estimated 
flight path angle 
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The circular arc angle 9 is. 


^ = 9 + i + X 


9 


'I' 



- X 


Now, we need to compute a command flight path angle, X^(0, V, U) , which 
will result in a circular (inertial) trajectory. Consider Fig. 38. Let 

the current position. O', be the origin of an inertial coordinate 
system. The problem is then to find x, y corresponding to the inter- 
section of Vg and loci, referenced to the known flight path tangent 

* 

through the current position. 

The wind vector anchor point locus is 
X = -y tan 9 - V 
and the aircraft velocity locus is, 

x2 = u2 - y2 

Combining: 

x^ = y^ tan^ 9^4* 2 Vy tan 9 + - y^ 

sec^ 9 y -h 2V tan 9 y - (u^ - V^) = 0 
and the solution to this is: 


2V tan 9 i /4 V2 tan2 0 -f 4 sec 2 0 (u2 - v2) 
2 sec^ 9 

V tan 9 + /tt2 9 - 

sec^ 9 


*Note: In general, there will be two solutions (clockwise and 

counterclockwise); therefore, suitable logic must be employed later. 
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Figure 38. Command Flight Path Angle 
for Circular Path 



cos 0 


(v sin 9 



and 


= -y tan 9 - V 

By pointing the aircraft (located at O') toward x^> the resultant 

velocity will be along the tangent, exactly. 

The command flight path angle corresponding to this is 


A a - 

= cos ^ X t 

llVall 

where t is the unit tangent: 

t = - sin 9 i + cos 9 j (for clockwise flight) • 
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Therefore, 


(- sin 0 + cos 0) 

= cos”^ _ _ - 

Jl . 2 

/xc + Yc 

We need to compute the command roll angle which corresponds to this path 
angle. To do this, first consider the ground speed along the circular 
arc and, knowing this, we can compute the bank angle required for circu- 
lar motion at radius R. 


'gc 


+ V, 


w 


= [- U sin (0 + X^) + Vjl 


+ U cos (0 + X^)J 


Therefore, 


ii'ec" ■ V 


= /u2 Sin2(0 -f X^) - 2 U V sin (0 + X^) + cos^ (9 + X) 


Ug^ = /u2 - 2 U V sin (9 + X^) + 


Now, consider Fig. 39. To generate the required acceleration 
2 

(* Ug^/R) directed toward the center of the circle, both longitudinal 
and lateral components must be generated, where 


U 


gc . > 

— sin X^, 
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Figure 39. Resultant Acceleration Vectors 

' ""<= X 

ay = — cos 


The means for generating accelerations in the aircraft axes are engine 
thrust and coordinated roll-yaw. Therefore 


= ma. 


= m 


8c 


sin X. 


tan 


-1 


U 


— cos X^ 
gR ^ 


If only roll commands are used, the situation is more difficult. In 
general, since the centripetal acceleration vector is skewed to the 
radius during crosswind flight, the center of curvature will translate 
relative to the desired center. The translation direction will be 
crosswind (W ^ % direction), with a magnitude bounded by (2R sin X^) 
per orbit (this is considerably less than if no wind correction was 
used) • 

To compensate for this, we could modulate the roll command on the 
up- and down-wind portions to give more or less lateral acceleration 
than originally, so the net center translation would be zero* This 
probably would give an elliptical, fixed orbit. 
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^ that 


To do this we would need to find a <(>^ =» 


^2it ^ 

I 

•'n 


d0 = ~ tan“l 

gR 


flu 

J 

•'n 


sin cos d9 


subject to the appropriate boundary conditions. This is an inconvenient 
task at best and a good guess might be: 

,2 




, "Sc 

^ sin X^ cos X^ sin 9 


(Approximate, for roll control only) 


This gives no change in acceleration at the crosswind positions (0 = 0, 
7t); an acceleration decrement at upwind (0 = tt/ 2), and increment at 
downwind (0 = 3 tt/ 2), of the proper magnitude. 

3 * Vehicle Motion 

A simplified set of fixed airspeed kinematic equations can be used 
to compute the resulting trajectory (for preliminary validation), i.e., 

X = U cos 

y = U sin ^ 

t S 

V = y tan (|)^, 


4* Results 

Figures E-1 to E-5 of Appendix E show resulting orbiting time histo- 
ries for the circular algorithm. These give x, y inertial positions 
(XI, X2) and heading angle (X3) as a function of time, for various vehi- 
cle airspeeds for a fixed 100 ft/sec wind velocity. As can be seen, 
when airspeed is high relative to wind speed, the time history is 
stable. As airspeed becomes closer to wind speed, the x and y peaks 
begin to drift away, linearly with time. Also, sharper "cusps” appear 
in the y trace, corresponding to higher turn rates on the downwind legs. 
When the wind speed is half the airspeed, the algorithm diverges, and 
orbit is not achievable. 
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These results are summarized In Fig. 40. This shows the xy tra- 
jectory as a function of wind-to-airspeed ratio. When the wind is less 
than 0.25 of the airspeed, the drift during the first orbit is quite 
small. Above this ratio the trajectory becomes increasingly elliptical 
and eventually diverges. The progression in drift rate suggests an 
exponential divergence with wind to airspeed ratio. 



Figure 40, First Orbit Trajectories for Circular Algorithm, 
as a Function- of Wind to Airspeed Ratio 

This is confirmed in Fig. 41 which shows the orbit drift velocity as 
a function of wind speed (both normalized by airspeed). In a realistic 
worst case situation where the wind aloft was 0.40 of the airspeed, the 
drift velocity would be one order of magnitude less than the wind speed. 
For a 15 minute flight of a low speed aircraft (V^ = 200 ft/sec) , this 
could still give nearly a 1-1/2 mile drift (on an extremely windy day). 
On a more typical day (V^^ = 30 ft/sec) with an approach configured 
HiMAT-llke vehicle (V = 300 ft /sec), the total drift would be less than 
50 ft! So this algorithm appears extremely sensitive to relative wind 
and air speed. 

The divergence of this circular algorithm is of course a consequence 
of the residual acceleration during portions of the orbit. As 
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mentioned, either incremental roll angle command modulation or auxiliary 
throttle control could be used to cure this. 

Further simulation runs indicated that this algorithm is insensitive 
to the initial orientation relative to the wind vector, and this factor 
is not related to the divergence. 

C. DEAD RECKONED WING-POINTING ALGORITHM 

A simpler and more effective way of maintaining a stable orbit can 
be realized by taking a few hints from these results as well as extant 
naval search and rescue flying techniques. Such "wing pointing** tech- 
niques give proper alignment of the centripetal acceleration vector, and 
are usually combined with altitude or airspeed modulation to give the 
desired stable orbit. 

In the current problem, altitude or speed modulation is not desir- 
able. Also, of course, a circular resultant orbit is an overly restric- 
tive requirement. Too, the divergence of the above orbiting algorithm 
may be related to failure to establish an imaginary orbit reference 
point. Consideration of all these factors leads to a simple, **dead 

reckoned wing-pointing** algorithm. 

1. Problem 

Fly closed path over ground, under steady wind, such that the 
acceleration vector lies along the wing, is of appropriate magnitude, 
and is directed toward a fixed point. 

Assume: 

• Constant airspeed magnitude, [V^l 

• Zero sideslip, no vehicle dynamics 

• Heading angle, 4^, available 

® Orbit begins at t=0"^ 

• Constant wind vector, 

• Fixed point on ground is x^, y^ 
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2* Analysis 


Vw 



yr 


» Vr 


Figure 42. Dead Reckoned Reference Form 

Let X = V =0. Dead reckon X , Y from 
o o r r 

Xj. = U cos i}; + 

= U sin ^ 

Then, the orbital reference angle is: 



From the definition of wing pointing, set: 


Then, 


^ = 


d^ 

dt 


a 





(33a) 

(33b) 


1 + u 


du 
2 dt 


, ^r 

where u = — 


iu 

dt 


(for + 7 t/2) 

_d ^ (xrYr - Yr^r) 

dt ~ 2 
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dll' 

■3t 


(xr7r - yA) 


[1 + (y^/x^)^] 


Since, 


xryr - yr^r 
(xr + yl) 


Xj. = U cos ij; H- 


w 


Then 


yj. = U sin ^ 


dt 


U sin 4^ y^. (u cos 4 ^ + V^) 


(x^ + yr) 

which leads to the command roll angle. 


(34) 


g tan (j)(, = iJi^U 

‘t'c = g 4’c (35) 

Figures F-1 to F-3 of Appendix F show the resulting time histories 
for this algorithm. This confirms that orbit stability is Independent 
of relative air and wind speeds (up to the limiting case of * 

The drift rate is zero, within the limits of numerical computation 
accuracy* The shape of the orbit itself is dependent on relative 
speeds, and results in a superellipse for non-zero windspeed* This is 
elongated along an axis perpendicular to the wind, as shown in Fig. 43. 
As wind speed increases, the dimensions of the superellipse expand. 
Also shown in Fig. 43 are the arc lengths covered by the vehicle during 
equal Increments of time. 
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0 



Figure 43. Effect of Wind and Airspeed on Vehicle 
Orbit Shape, Wing Pointing Algorithm 


Further analyses have shown that this algorithm obeys Kepler's Law 
of Areas, i.e., equal areas of arc are swept out in equal lengths of 
time. Also, it is not an inverse distance-squared form resulting in 
simple planetary-like ellipses, but instead results in a higher order 
trajectory. 

D. CONCLUSIONS AND RECOMMENDATIONS 

A wing pointing algorithm, where the heading rate is varied so that 
the wing always points towards a dead reckoned orbit reference point, is 
the recommended approach to maintaining autonomous RPV orbits in the 
presence of winds. The resulting superelliptical orbit is entered when 
the vehicle is parallel to the estimated wind. The required heading 
rate is translated to a roll angle command. The basic requirements for 
mechanizing the algorithm are: onboard inertial heading and airspeed 
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signals; an estimate of wind speed at orbit initiation; and onboard 
firmware Aiardware computation of the needed trigonometric, arithmetic 
and integration operations (of Eqs. 33, 34, and 35), 

Clearly, the next step in these analyses should be to examine the 
effects of non-idealized assumptions, such as: 

• Actual wind variation (magnitude and direction) 

• Sensor noise, drift and resolution 

• Vehicle dynamics 

• Interface with other autonomous modes 

More real-world data is needed to perform these analyses. 
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SECTION VI 


AUTONOMOUS NAVIGATION AMD LANDING CONSIDERATIONS 


Upon loss of uplink signal it is desired (Ref. 3) that the onboard 
control system autonomously navigate from any point at which uplink loss 
occurs to a safe orbit area, deorbit, and land. The performance 
requirements for the system include maintenance of a position accuracy 
within a CEP of 100 feet when operating with radar updates. For safe 
orbiting, the required accuracy is a CEP of 3 to 4 nmi. An additional 
performance objective is a CEP of 1 nmi as this would allow considera- 
tion of an autonomous recovery mode (landing) as an extension to the 
system concept. The system is to operate in a latitude/longitude coord- 
inate system with radar position updates in the same form. 

This can be accomplished by using a simple navigation and control 
scheme similar to that employed in the Aqulla RPV (Ref. 14). The Aquila 
is a small reconnaissance drone which includes an autonomous dead reck- 
oning and orbit mode upon loss of telemetered uplink command. 

A key element in the Aquila system is a strapdown magnetometer which 
senses vehicle heading. The magnetometer and associated electronics is 
off-the-shelf, small, and light weight (4.5 cubic inch, 3.5 oz.). It 
measures the magnetic field in vehicle x, y, z coordinates with the x 
and y components used to calculate magnetic heading. When operating in 
a relatively small area (such as the DFRF test area) of known declina- 
tion and magnetic anomaly and with calibration for vehicle magnetic 
fields, the sensed heading is quite accurate. 

Under normal operations, vehicle range, azimuth, and altitude from 
the FPS— 16 tracking radar can be uplinked and stored in an onboard navi- 
gation computer module. Vehicle position and groundspeed can then be 
computed in either polar or rectangular coordinates. When combined with 
onboard sensing of airspeed and heading, the direction and magnitude of 
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wind can be determined and maintained in the computer. Likewise » head- 
ing and range from the vehicle current position to the position coordi- 
nates of a safe orbiting area can be computed continuously. Then upon 
loss of uplink the vehicle can be immediately commanded to the heading 
of the safe orbit area. Estimated groundspeed can be calculated from 
measured airspeed and the wind data stored in the computer. The vehicle 
is then dead reckoned to the safe orbit position on the basis of head- 
ing, groundspeed, and elapsed time. Again, the basics of this naviga- 
tion system are contained in the Aquila RPV and need only minor computa- 
tional additions. 

Upon reaching the coordinates for entering the autonomous orbit 
mode, the computer can command the vehicle into the wind pointing orbit 
of Section V and a barometric rate of descent. The vehicle will then 
spiral down until a deorbit command can be given. Deorbit is estab- 
lished by radar altimeter and heading measurements. After reaching a 
preselected radar altitude and the appropriate heading, the vehicle is 
commanded to a wings level, constant heading descent to the lakebed run- 
way. Autonomous flare and landing can then be accomplished by the pro- 
posed (or present) flare control law (Section III). 

Based upon available data, it appears that the Edwards precision 
bombing range would provide a safe orbit area with the possibility of an 
autonomous landing. Figure 37 is a map of the precision bombing range 
and shows that an orbit of up to 5 miles diameter can be easily accommo- 
dated. The orbit can be made target to an extension of the runway 
25 centerline. By flying the orbit in a clockwise direction the vehicle 
can intercept the runway alignment and then fly a straight approach to 
landing. 

Figure 38 is a topographical representation of the terrain elevation 
versus range along the approach to runway 25. It shows the ground to be 
essentially a constant slope of approximately 0.7 deg. If the Fig. 2 
2.8 deg descent is retained, the approach would be flown at about 2 deg 
with respect to the ground. 
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Figure 37, Proposed Autonomous Orbit and Approach Pattern Within the 
Precision Impact Range Boundaries 
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Figure 38. Ground Elevation on Runway 25 Centerline Approach 
Through Free ision Impact Range 



This approach may or may not result in a landing on runway 25 but, 
in any case it should 

— be close 

— not cause any safety problems 

— minimize risk of property or vehicle damage 

Thus, it appears that the autonomous navigation, orbit, and landing 
is not only feasible but highly promising. 
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APPENDIX A 


A. RECONSTROCTION OF THE PA-30 LONGITUDINAL AIRFRAME 

AND FLIGHT CONTROL SYSTEM M3DEL 

This appendix documents an effort to reconstruct the mathematical 
model of the PA-30 airframe and flight control system for the longitudi- 
nal modes employed in Ref. A-i. There Edwards presents flight test data 
for evaluation of the HIMAT approach and landing backup control system 
(BCS) using the PA-30 airplane, results of a nonlinear simulation imple- 
mented on the NASA Dryden CYBER computer, and results from linear analy- 
ses of the pitch autoland system. It was desired to use these analyti- 
cal models in this study (Ref. A-2). However, a check of the Ref. A-1 
analysis showed some large discrepancies. It was then necessary to 
reconstruct the Ref. A-1 model and this required combining the Ref. A-1 
data with data from other sources, i.e., Ref. A-3 through A-7. As a 
result, a linear, small perturbation longitudinal model of the airframe 
and BCS control system has been developed which matches selected root 
locus and Bode plots from Ref. A-1. This model is presented in Subsec- 
tion 1 and the process by which it was generated is discussed in Subsec- 
tion 2. 

1. Reconstructed Model 

The linear longitudinal model of the airframe and flight control 
system presented in Ref. A-1 is reproduced here as Table A-1. The pres- 
sure altitude of 3350 ft above sea level corresponds to a radar altitude 
of 1000 ft above ground level. The true airspeed of 186 ft/sec is an 
approximation of the approach speed of 110 kt indicated. The values of 
and N^ imply a flight path angle of -5 deg and at this speed is con- 
sistent with a rate of descent of 1,000 ft/min. These flight parame- 
ters, along with geometric and inertial data obtained from Refs. A-3 
through A-7, are summarized in Table A-2. 
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Linearized small perturbation longitudinal airframe equations are 
given in Table A-3 consistent with the assumptions and notation of 
Ref* A-8 for stability axes* The nondimensional and dimensional stabil- 
ity and control derivatives finally evolved are given in Table A-4. 
These derivatives when used with the equations of motion of Table A-3 
produce the Table A-5 bare airframe transfer functions* The Table A-5 
denominator checks closely with that of Table A-1* Table A-1 does not 
contain transfer function numerators (zeros) but the correctness of the 
Table A-5 zeros are implied from the loop closures to be discussed 
below* 

Figure A-1 shows the portion of the BCS flight control system con- 
sisting of feedbacks to the elevator. (The autothrottle system will be 
considered later). The control system gains are summarized in Table A-6 
as a function of switch position* It should be noted that the gains in 
the proportional and Integral pitch rate feedback paths are not the same 
as given in Ref* A-1. Transfer function models for key control system 
elements are presented in Table A-7. Provision is made in Fig, A-1 for 
effective time delay and hysteresis in the elevator feedback but they 
were not included in the initial analysis* 

To verify that the model outlined above matches the closed-loop 
dynamics of Ref* A-1, loop closures were made in the same sequence as in 
Ref. A-1* Figure A-2 shows • the closure of the rate of climb loop, 
h-^6e> with the a^ and q loops open. Without the latter loops closed, 
the phugoid mode is driven unstable* Figure A-3 s];iows the comparable 
locus of Ref. A-1. The close agreement verifies in Table A-5 and 
the use of KIHD in Table A-6* 

Figure A-4 shows the closure of the pitch rate loops about the rate 

of climb loop but with the a^ loop open* The pitch rate closure resta- 

bilizes the phugoid mode while destabilizing the servo mode. This 

closure corresponds to the loci labeled "q^^” in Fig. A-5 (taken from 

Ref. A— 1). The close agreement of Figs. A-4 and A-5 at switch position 
0 

4 verifies the numerator and the KQAL and KIQ gains shown in 
Table A-6. (As will be explained in Subsection 2, these q loop gains 
are not those given in Ref. A-1 but rather come from Ref* A-9). 
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Finally, Fig. A-6 shows the closure of the last loop, 
the other loops closed* It may be seen that this loop has rather mini- 
mal effect. This closure should be compared to the 
Fig. A-5* The close match of the final root location requires the B gain 
set at 4 (KINZ = 0.1919 deg/g) rather than at 2 as in Ref. A-1 . Since 

this loop has little influence this difference can be neglected. The 

a^ 

results tend to verify the zeros of Table A-5. 

The Fig. A-7 solid line Bode plot shows the resulting open-loop 
transfer function with the loop opened at the servo command point (the 
dashed curves will be discussed in the next section). The system gains 

correspond to the 441 setting for Table A-6 and thus this response may 

be compared to Fig. 24 in Ref. A-1, reproduced here as Fig. A-8. This 
comparison shows a good match between the amplitude and phase curves at 
low frequencies. The comparison does not appear good at high frequency 
however, because Fig. A-8 reflects a discrete frequency response in the 
w^-domain and includes digital artifacts not included in the s-domain 

analysis of Fig. A-7. As noted in Fig. A-7 this phase lag difference is 
equivalent to an effective time delay of about 0.063 sec which is appro- 
priate for the uplink delay of Fig* A-1. 

2. Reconstruction of the Model 

Because complete details were not provided in Ref. A-1, considerable 
iteration and trial and error gain analysis was required to obtain the 
match presented in Subsection 1. The altitude and speed were specified 
for the reference case as 1000 ft AGL and 110 kt indicated. The flight 
path angle, Y> was derived from the BCS scheduled altitude rate command 
of 1000 ft/mln at a radar altitude of 1000 ft. This is consistent with 
the n and n values given in Table A-1 since, if we assume the use of 
stability axes, 

sin~^ (0. 0849) = cos“^ (0. 9964) = 4.9 deg 


TR-1 180-1 


A-3 



From this flight condition the trim may be estimated using 


L = W (A-la) 

T = D + WYo (A-lb) 

which gives an estimate of the trim lift coefficient as 


Cl = 


w 

(l/2)pUoS 


= 0.491 


(A-2) 


Since airframe aerodynamic, mass and geometric data are not specifi- 
cally identified in Ref. A-1, an initial model was constructed from data 
in Refs. A-3 through A-7 and used to generate characteristic roots which 
were then checked against those shown in Table A-1. This initial esti- 
mate resulted in real short-period roots considerably different from 
those of Table A-1. After exploring the possibilities for uncertainties 
in the model, it was found necessary to extract the airframe parameters 
from the matrix equations of Table A-1. The state vector implies that 
the first four equations of the 8x8 "A" matrix are the longitudinal 
equations. It should be noted that this matrix indicates coupling of 
sideslip into the longitudinal motions and coupling of angle-of-attack 
and speed perturbations into the lateral-directional dynamics. The 
coupling of sideslip into the longitudinal motion does not seem justi- 
fied since the airplane is supposedly in a symmetric (8q “ 0) flight 
condition. The fourth equation may be seen to be 0 = q consistent with 
Table A-3. The third equation is the pitching moment or q equation as 
indicated by the fact that the attitude element, is zero. The first 
equation is the x-force equation since the 0 element is consistent with 
(-g cos Sq)* The fact that the q term in this equation is zero indicates 
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the use of stability axes# Finally, the second or a equation is equiva- 
lent to the z-force equation of Table A-3 divided by as indicated by 
the fact that the q term in the second equation is exactly 1. Thus, the 
equations used in Ref. A-1 can be considered equivalent to those of 
Table A-3 with the exception of the state variable formulation and the 
use of a instead of w. 


It may be immediately noted that the Ref. A-1 formulation neglects 
the derivative since its presence would require inclusion as an 
apparent product of inertia term in the inertia matrix. This apparently 
was not done in Ref. A-1. Neglecting reduces considerably the 
modeled airframe-alone short period damping. 


Since the equations of motion are divided through by mass and moment 

of inertia, it is not possible to verify these parameters independently 

of the aerodynamics. However, assuming the values obtained from 

Refs. A-3 through A-7 are correct, the aerodynamic stability derivatives 

may be identified in the Ref. A-1 matrix equation elements. From the q 

element of the pitching moment equation, = -3.38 sec ^ which 

implies = -14.4 rad"^, considerably lower than the value of 

-25 rad'"'*' given in Refs. A-3 and A-5. The a element of the a equation 

is “ 1.98 sec"^ and indicates a Cl^ of approximately 6.1 rad 

This is about 26 percent higher than the value obtained from the wind 

tunnel data of Ref. A-6 ^and probably reflects the 30 percent increase in 

lift curve slope which was mentioned in Ref. A-1 as being required for 

the PA-30 simulation model. The M derivative as extracted from the a 

a 

element of the pitching moment equation is —5.68 rad which implies 
Cm “ -0.325 rad“l , and, when used with the value above, indicates a 
static margin of only 5 percent c, which would appear to be too low for 
actual flight of the PA-30. Some additional differences were also found 
in the drag and elevator effectiveness derivatives. 


Once the airframe parameters were extracted, loop closure checks 
were performed as discussed in Subsection 1. Originally these root loci 
were generated using a gain/gain switch table in Ref. A-1. However, the 
Ref. A-1 A gains (q loop gains) did not produce a good match with the 
Ref. A-1 loop closures. It was found that q loop gains given in 
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Ref* A-9 provided an excellent match and these are given in Table A-6 
combined with B and C gains from Ref. A-1. It is further assumed that 
these later of loop gains were actually used in flight. 

A comparison between the derivatives extracted from the Ref. A-1 
matrix and those derived from the Refs, A-3 through A-7 wind tunnel 
data, theoretical estimates, and simulation models are compared in 
Table A— 8, The resulting difference in system dynamics is shown in 
Fig. A-7 where the frequency response shown by the solid line corre- 
sponds to the reconstructed model as derived from the Ref. A-1 matrix. 
The amplitude and phase curves shown in the dashed line corresponds to 
the Table A-8 derivative set obtained from Refs. A-3 through A-7. While 
the general shape of the magnitude and phase curves is quite similar, 
the Refs. A-3 through A-7 derivatives produce a much lower gain over the 
frequency range up to 10 rad/sec. This would result in considerable 
droop and sloppy control in the region above 1 rad/sec. Since the 
Ref. A-1 simulation appeared to match flight results very well it is 
presumed the Ref. A-1 aero coefficients are the better set. Thus, while 
the stability derivatives reconstructed from the Ref. A-1 matrix will be 
used in further analyses, it should be remembered that there is signifi- 
cant sensitivity to uncertainties in the airframe dynamics that may 
effect closed-loop response. Furthermore, it is not obvious that the 
Refs, A-1, A-3, or A.— 5 simulations adequately model the propwash 
effects. 

3. Airframe Characteristics at 80 kts 

After the airframe characteristics were reconstructed for 110 kts, 
the stability derivatives were adjusted for n^ = 1 g, y = 0 flight at 
80 kts (see Table A-9). The four elevator transfer functions for this 
case are tabulated in Table A-10. 
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TABLE A-1. PA-30 APPROACH AND LANDING MATH MODEL FROM 

REFERENCE A-1 (REPRODUCTION OF TABLE 2 IN REFERENCE A-1) 


H * 

3350«0 

V « 185.00 

iACH = 

.17 

THRUST a 383. 

124 


ALP a 

4.8697 

8TA a 

0.Q030 

THA * 

4. 8697 

PHI * 

0.0000 

OE » 

.6169 

OA a 

0.0030 

OR a 

0. 0000 

00 a 

Q.OQQQ 

P ■ 

O.QOC8 

Q = 

0.0000 

R a 

0. 0000 



NX * 

.0 849 

NY a 

-.0646 

NZ * 

-.9964 




V. ALP.Qf THA 
V 

,PtR.8TA.PHI 

a 

q 

e ' p 

r 

y9 

4 > 

-.3570E-31 

.1492E»02 

0. 

-.3216E»02 [ 0. 

0. 

-«2204E»00 

0. 

-.1777E-C2 

-.1978E»01 

« lOOOEf 01 

• 3309E-12 1 0. 

0. 

C . 

0. 

-.376ie-03 

-.5679E»01 

-.3378Ef01 

0. (0. 

0. 

-.4446E-03 

0. 

C • 

0. 

.1C00E»*01 

0. Jo. 

0. 

0. 

0. 

0. 

-.73l3E»00 

C. 

0. ] -.4105Ef01 

. izieeo'i 

-.40056^01 

0. 

0. 

-.1262Ef0Q 

0. 

0. I -.5895E^00 

-. 5033Ef 00 

.37C7Ef0l 

0. 

.6441E-05 

-.7697E-&2 

c. 

0. [ .84R9E-01 

9964€f DO 

-.1029EfOO 

.1732E^00 

0* 

0. 

0. 

0. I .lOOOEfOl 

. 8520E-Q1 

0. 

0. 

OEtTHRUSTtDAtCR 

Se T 

Sa 

Sr 




.5q55£>31 

.8618E-02 

0. 

0. 




-.2618Ef00 

-.2470E-04 

0. 

0. 




-•4020£»02 

♦9080E-04 

0. 

0. 




0. 

0. 

0. 

0. 




0. 

-.1955E-03 

•66246^01 

.7575E^00 




0. 

-.1941E-03 

.1670EfGQ 

-.3028Ef01 




0 • 

-.3110E-05 

0. 

.20Q9E-01 




0. 

0. 

0. 

0. 




NZtNY 







-.1028E-01 

-.1129E»02 

0. 

•1639E-04 0. 

0. 

0. 

0. 

•1388E-13 

-.4426E-01 

0. 

0. . 0. 

0. 

-.5916E»00 

-.1924E-03 

-. 1485E»01 

-.1188E-03 

0. 

0. 




0. 

-.1788E-04 

0. 

. 1155Ef00 





MACH 

NO. a 

.17 

ALTITUDE a 

3350.00 ET 


DENOMINATOR 

ROOTS 





REAL 


IMAGINARY 

OMEGA 

ZETA 

TAU 

-2.67980 

-2.67980 

-.01591 

-.01591 


-2.28068 
2.26068 
-. 155C2 
.15502 

3.51893 
3. 51893 
. 15583 
.15583 

• 76154 
.76154 
. 10213 

• 10213 

0.00000 
0. 00000 
0. QOOOO 
0. 00000 
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TABLE A-2. FLIGHT CONDITION AND AIRCRAFT 
CONFIGURATION 


Altitude, h 
Airspeed, U^ 


Flight path angle, Yq 
T rim, C^ 

Trim, Cj^ 

Reference planform 
area, S 

Reference chord, c 

Weight, W 

Pitch axis moment 
of inertia, 


1000 ft above ground level 
3350 ft above sea level 

186 fps true (approximation of 
110 KIAS) 

-5 deg 

0.49 

0.056 

178 ft^ 

5.00 ft 
3600 lb 
1900 slug-ft^ 


TR-1 180-1 


A-8 


TABLE A-3. LONGITUDINAL AIRFRAME EQUATIONS 


• Linearized small perturbation form 

• Stability axis, Yq = , Og 

Equations of Motion 

(s - Xg)u - + g cos GqQ 

— ZgU + (s ~ + (~UqS + g sin 0o)9 

-MgU - (M^ + Mj,)w + (s^ + Mqs)0 


0 


^6e<Se + X5q,6x 
M6e«e + 


S0 = q 


Altitude Rate at the h 

sh = h = 


sin 9 qU - cos 0 qW + Uq cos 


z-Axis Acceleration at Accelerometer Location* a^ 

a 


a 2 = w - Uq 0 + g sin 0g0 - ^ 0 

cL 


Accelerometer is assumed to be mounted at the instant center for 

elevator inputs, = 1,21 ft forward of c.g. 

a 
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TABLE A-4. STABILITY AND CONTROL DERIVATIVES (110 kts) 


NON- 

-DIMENSIONAL 

DIMENSIONAL 



STABILITY 

DERIVATIVES 


Cd 

*^u 

0 

Xu 

= 

t-«D - C„„) 

-0.03567 sec“^ 

Cd 

0.24 rad“^ 

Xw 


Ccl - OdJ 

0.07962 sec"^ 


0 

Zu 

■ = 

f (->=L - Cl^) 

-0.3121 sec~l 


6.15 rad"^ 

Zw 

= 

- ‘=») 

-1.976 sec“^ 

Cm 

0 

Mu 

= 

pSUc , , 

ly (Cm + Cm^) 

0.0002425*(ft/sec)"^ 


i 

-0.325 rad~^ | 


= 

pSUc 

2Iy ^M„ 

-0.03048 (ft/sec)”^ 

Cm^ 

0 1 

MvJ 

* 

pSc^ 

0 ft~l 

Cm 

”q 

-14.4 rad"^ 

Mq 

= 

pSUc^ 

«y '’"q 

-3.376 sec“l 

CONTROL (ELEVATOR) DERIVATIVES 


0 rad ^ 

X6e 

= 

pSU^ r 

- ^ 

0 ft/rad sec^ 


0.818 rad“^ 

"«e 

= 

- Ct 

2m 

-48.45 ft/rad sec^ 


-2.30 rad"^ 

Mde 

= 

pSu2 

21 

Oq 

-40.12 (rad/sec^) ^ 


it 

Due to thrust axis 0.75 ft above and parallel to 
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TABLE A-5. AIRFRAME -ALONE ELEVATOR TRANSFER FUNCTIONS (110 kts) 

A = [0.1185, 0.1620] [0.762, 3.51] 

N6 = -3.86(3.59)(-179.8) ^ 

e rad 

= -48. 5(157. 4)[0. 1086, 0.231] 

= -40.1(0.0488)(1.926) ^ 

rad 

N^^a = 0.0864[0.661, 0.0271] [-0.982, 409.] 

n5 = 48.3(0.0356)(-15.64)(19.03) ^ 

e rad 

Notation: K(a)[c>u)] = K[s + a)[s^ + 2501 s + u>^] 
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TABLE A-6. CONTROL SYSTEM GAINS 


SWITCH 

POSITION 

GAIN 

A 

B 

c 

KQAL 

rad/ (rad/ sec) 

KIQ* 
rad/ rad 

KINZ 

deg/g 

KIHD 

deg/fps 

1 

0.1 

0.3696 

0.0768 

0.275 

2 

0.15 

0.5544 

0.1151 

0.275 

3 

1 

0.2 

0.7392 

0.1535 

0.275 

4 

0.25 

0.9240 

0.1919 

0.275 


* 

KIQ = 3.696 KQAL 


The gains for each of the 3 loops (q-»-6e» be varied 

independently by setting the "A", "B", and "C" switches. The gain set is 
specified by the ABC switch sequences, e.g., the "441" gains are (KQAL = 
0.25, KIQ = 0.9240, KINZ = 0.1919, KIND = 0.275) 
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TABLE A-7. CONTROL SYSTEM ELEMENTS 


Datametrics Rate-of-Climb Sensor 


Gh(s) 



Elevator Servo 



1 




+ 


2(0.7) 
26.6 ® 


+ 1 ] 


No Delay or hysteresis 



1.0 
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TABLE A-8. COMPARISON OF STABILITY AND CONTROL DERIVATIVES FROM 
REF. A-1 MATRIX AND FROM REFS. A-3 THROUGH A-7 (110 kts) 


DERIVATIVE 

FROM REF. A-l* 
MATRIX 

FROM REFS. A-3 
THROUGH A-7 
(Wind Tunnel Data, 
Theoretical Estimates) 

^Lcx 

1/rad 

- ... . . - 

6.15 

4.87 


0.056 

0.047 


1/rad 

0.24 

0.334 


1/rad 

-0.325 

-0.974 


l/rad 

0. 

-14.6 

% 

1/rad 

-14.4 

-25.0 


l/rad 

0.818 

1.03 


l/rad 

-2.30 

-2.64 


Reference A-1 matrix 
here are the same as Table 


reproduced in Table 
A-4. 


A-l, 


derivatives shown 


TR-1 180-1 


A-l 4 





TABLE A-9. STABILITY AND CONTROL DERIVATIVES (80 kts) 



NON- 

-DIMENSIONAL 

DIMENSIONAL 




STABILITY 

DERIVATIVES 




0 

Xu 

= 

(-<^D - C„J 

-0.0381 sec"l 



0.46 rad"^ 

Xw 

= 

(-=1. - <=D„) 

0.112 sec"l 



0 

Zu 

= 

'’f (-<=L - C.^) 

“0.443 sec""^ 



6.15 rad ^ 


= 

(-^L„ - C„) 

“1.48 sec“^ 



0 

Mu 


fc + C 1 

0.00168*(ft/sec)"l 






-0.325 rad"^ 

My, 

= 

pSUc 

2Iy 

-0.0228 (ft/sec)~l 


Cms 

0 

Mw 

= 

pSc^ 

4ly 

0 ft“l 


Cm 

Mq 

-14.4 rad"^ 


= 

pSUc^ 

«y % 

-2.52 sec~^ 


CONTROL (ELEVATOR) DERIVATIVES 

-- 


0 rad~^ 

X6e 

= 

- c.^ 

2m 

0 ft/rad sec^ 



0.818 rad"^ 

"^e 

= 

PSU2 

2m 

“26.3 ft/rad sec^ 

- 


“2.30 rad"^ 

M6g 

= 

PSU2 

2Iy ^Mg^ 

“21.8 (rad/sec^) ^ 


Due 


to thrust axis 0.75 ft above and parallel to 
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TABLE A-10. AIRFRAME -ALONE ELEVATOR TRANSFER FUNCTIONS (80 kts) 

A = [0.0528, 0.243] [0.767, 2.62] 

N5 = -2.94(2.74)(-126.6) ^ 

e rad 

Nfi = -26. 3(114. 3)[0. 0577, 0.322] 

e rad 

Nfi = -21.8(0.0741)(1.420) 

°e rad 

Ns^ = -26.3(0)(-0.000583)(-11.54)(14.10) 

rad 

• 

N* = 26.3(-0.000583)(-11.5)(14.10) ^ 

e rad 

Notation: K(a)[i;,uj] = K[s + a)[s^ + 2cos + uy^] 
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HDLCMS 



A-17 






Figure A-2. Closure of the h->-6e Loop With a^, q'*'<Se 
and u+iS']’ Loops Open (110 kts) 
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ORIG^AL PAGE IS 
OF. QUALITY 



Figure A— 3 . h Root Locus With Datametrics Pressure Sensor, 
3,350 ft, 110 kts (Reproduction of Figure 18, Reference 1) 
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page is 
op POOR QUAtlTV 


Servo 



Gain Switch = I 


J 






Figure A-5. q and Root Locus With Datametrics Pressure Sensor, 3,350 ft 
110 kts, c Gains, Run 23 (Reproduction of Figure 21, Reference 1) 
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APPENDIX B 


PA30/BCS ADTO THROTTLE ANALYSIS 


!• Auto Throttle Mechanization for the PA--30 System 

A complete analysis of the existing PA-30 auto throttle mechaniza- 
tion could not be made because of uncertainties regarding throttle 
deflection units and the throttle effectiveness derivative, X 5 X. This 
situation produced an uncertainty in the overall loop gain for the auto 
throttle* Further, there was no linear auto throttle analysis in 
Ref. B-1, such as for the elevator loops, which could be used to resolve 
these questions. However, some approximate analysis has been done to 
understand the gross effects and to approximate auto throttle behavior 
sufficiently for use in analysis of the more important elevator loops. 
The auto throttle as shown in Fig. B-1 is basically an airspeed command 
system consisting of a U+ 6 x feedback and an airspeed command, U^, 
scheduled with radar altitude. The airspeed schedule presently used is 


! 74.4 + 0.0375 h^ kts 
(125.7 + 0.0634 fps) 
80.0 kts 
(135.2 fps) 


Decel region 


Flare 


(B-1) 


It may be seen that the auto throttle is a pure airspeed to throttle 
loop in the flare mode, but in the decel region there is in addition an 
effective radar altitude-to-throttle feedback. This additional feedback 
will be considered briefly to demonstrate that its effect on the air- 
speed dynamics is negligible and that it may be treated purely as a 
quasi-static airspeed command. 

For small perturbations about any (quasi-static) reference speed, 
the auto throttle dynamics may be linearized. If we consider the 
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Figure B-1. PA-30/BCS Auto Throttle 

closure of a throttle loop first before closure of the elevator loops, 
it may be shown (Ref* 8“2) that the throttle feedbacks modify the eleva* 
tor and gust transfer function denominators and numerators* 

2. Auto Throttle Effect on System Dynamics 

As indicated in Eq. B-2 


-1 

^ • u, h^‘^6'p 




(B-2) 


where 


X = u, 9, a, etc. 

<5 = iSg, Ug, Wg, etc. 
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From Eq. B-*l and Fig. B-1 


-0.0375 Kx 




Decel region 


6x unit /ft 


0.0 


Flare 


(B-3) 



-0.24.5s 
s I (s + 5)2 " 


1.2s + Kx(s + 5)2 
s(s + 5)2 


6-j; unit/kt 


(B-4) 


The above feedback transfer functions are written for negative feedback 
consistent with Eq. B-2. In addition to the uncertainty in the overall 
auto throttle loop gain, there is some uncertainty in the value of ICj, in 
Fig. B-1 (denoted as ”KIQCT" in the NASA DFRF documentation). In 
Ref. B-1 


67'unit 

Kx = "KIQGT" = 0.00218 0.00545 — (A-5a) 

In Ref. B-3 

6t unit 

Kx = ”KIQCT" = 0.00368 0.00920 — ■ ■ . — (A-5b) 


However, for purposes of approximate analyses, we may assume that Kj, is 
on the order of 10 . 
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3. Nature of the Numerator 

The numerator of the feedback transfer function in Eq, B~4 will 
be a quadratic with real roots* The nature of the roots, as a function 
of ICj. may be seen in the system survey sketched in Fig* B-2. This indi- 
cates that for K-p = 10*^, there is one large root greater than 10 rad/ 
sec and the small root which to an asymptotic approximation is 20.8 Kj, 
rad/sec* Thus, from Fig. B-2 and Eq. B-4 


Kp(l/Tf)(l/T^) 

(0)(5)2 

^ Kp( 20. 8Kp) ( 1 . 2/Kp) 6p unit 
(0)(5)2 


(B-6) 



a) Bode Roof Locus 


b) Conventional Root Locus 


Figure B-2. 


System Survey Sketch for 


the G^T 


Numerator 


B-4 


TR-1 180-1 


For the frequency region of interest which is well below the high fre- 
quency zero. 


^ 1.2(20.8 Kx) 6x unit 

■ (0)75)2 ~ 


which has an amplitude frequency response as sketched in Fig. B-3. It 
should be noted that the asymptote between the two zeros is independent 
of Kj,. 



log oj ( rad /sec) — ► 


Figure B-3. 


Frequency Response Sketch of 


Amplitude 


4. Relative Contribution of the Airspeed and Altitude 
Feedback Loops to A' in the Decel Region 

From Eq. B-2 it may be seen that the relative contribution of the 
altitude and airspeed loops may be compared from the ratio 


^uh 


4 ^ 


(B-8) 
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The contribution to the above ratio, which is dependent on the control 
system parameters is 


gJt -0.0375Kx(5)2 

6x ~ Kx(20.8Kx)(1.2/Kx) ft 

^11 


It may be seen from the asymptotic frequency response sketch in Fig* B-4 
that this ratio will be much less than 1 at all frequencies, i*e*, at 
least -28 db or less. The airframe numerator ratio in Eq. B-8 for the 
110 kt case is 


-0.000779(-0.0783)[0.0684, 13.62]/(0) ft 
u 0.00894/1. 69(0. 242)[0. 742, 3.44] kt 


(B-10) 



log (o (rad /sec) — ► 


Figure B-4, Frequency Response Sketch of 

The sketch of the Eq. B-10 frequency response shovm in Fig. B-5 indi- 
cates that the numerator ratio rolls off rapidly with increasing fre- 
quency and that the overall ratio given in Eq. B-8 will be on the order 
of 10"“1 or less in the region of interest for vehicle dynamics, thus 
indicating that the h-^6x feedback is significant only at very low fre- 
quencies, i.e., the steady state. Thus, as indicated earlier, the h^6'i^ 



effect can be treated purely as a quasi-static variation in the airspeed 
command. Numerical surveys of the various coupling numerator ratios 
indicate that in an analogous manner the effect of the h feedback may 
generally be neglected on the numerators of Eq. B— 2 — at least for the 
elevator transfer functions. 

5. Approximate Augmented Denominator 

From the preceding considerations, the effect of the auto throttle 
on the characteristic polynomial may be approximated as 


= A + 


(B-11) 
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There is some effect of the stability derivative (i.e., pitching 
moment due to throttle) which occurs as a consequence of thrust offset. 
This pitching moment contribution primarily effects the low frequency 
real root in the speed to throttle numerator. However, for an approxi- 
mate analysis of the auto throttle closure, both the and the 
effects may be neglected and thus 

^6^ “ ^i5js[s^ — + Mq + M^js + (Z^Mq — ] 


= X^^s Agp (B-12) 

Under these idealizations, the second order speed zero cancels the short 
period pole in the speed to throttle transfer function as seen in 
Eq. B-13 


N(5t ^ ^6'pS^sp 
^ ~ ^*^sp 

The phugoid mode is thus ideally the speed response mode. 

Thus, the closure of the auto throttle loop may be seen from the 
system survey of Eq. B-14 


X5q.s 

lCp,“p] 


(B-13) 


^u^ ^5^ . 1.2(20.8%) Xj^CO) 

^ “ (0)(5)2 * t4p,a)p] 


(B-14) 


shown in Fig. B— 6 as a Bode root locus. Figure B— 6 shows that the low 
frequency zero of the airspeed-to-throttle equalization (Fig. B-2) is 
positioned near the phugoid to produce a "K/s region" above the phugoid 
out to u) = 5 rad/sec. The auto throttle crossover frequency cannot be 
determined, of course, because of the uncertainty in loop gain noted 
previously. However, Fig. B-6 indicates that good auto throttle per- 
formance could be achieved with closures in the K/s region. Assuming 
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Figure B-6. Bode Root Locus Sketch of the Auto Throttle Closure 

that the auto throttle does have a reasonable gain, then it is reason- 
able to consider its effect as an ideal airspeed constraint in the speed 
response (phugoid) region. Assumption of a pure speed constraint leads, 
of course, to the short period model. Ref. B-2 for the airframe dynam- 
ics. Thus, even if the .auto -throttle bandwidth is lower than the short 
period mode frequency, the short period equations of motion should pro- 
vide an adequate auto throttle model for analyses of the elevator feed- 
back loops. 
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APPENDIX C 


PA-30/BCS ELEVATOR FEEDBACK LOOP ANALYSIS 

!• Short Period Model of the Airfrane Dynamics 

The auto throttle analysis, Appendix B indicates that even though 
the effective gain of the auto throttle (and hence the bandwidth) is 
unknown, it is reasonable to treat the auto throttle implicitly by use 
of the speed constrained, short period equations of motion for analysis 
of the elevator loop closures* The short period model, as given in 

Ref* C-1, pg. 637 is summarized below 


Equations of Motion: 

[s - Z^]sd + Zqj 9 = + Z^n 

(M^s + M^^)sd + [s^ - (Mq + M^)s - (C“l) 

Characteristic Function: 

A = s^{s^ “ (Mq + + Z^)s “ - Z^^Mq ) } = s^Agp (G-2) 

Numerator 

A ^5 

N6 = -Z5[s 2 - (Mq + m;)s - (m„ - — Z„) ] 

Nq = Zq[s^ ~ (Mq + M(j)s ~ ~ Z(j) ] 

N0 = s{(M 5 + ZfiM^js - (MsZ„ - ZfiMw)} 

Nq = -s { (Mf^ + ZqM^^js - (m^Z,, — } (C— 3) 
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2, BCS Path Control System — PA-30 


The BCS elevator feedback loop structure as presently Implemented on 
the PA— 30 is shown in Fig. A— 1 and consists of feedbacks of pitch rate 
from a pitch gyro, acceleration from an accelerometer, altitude rate 
from a barometric rate of climb sensor and altitude from a radar alti- 
meter. A linear analysis of this control system without speed con- 
straint, i.e., neglecting the auto throttle was given in Ref. C-2 and 
that analysis is reviewed in Appendix A. The elevator loops will be 
briefly analyzed in this section with the auto throttle effect included 
implicitly through the short period model. This analysis will provide 
some perspective for the proposed system to be discussed in Appendix D. 

3. Pitch Rate to Elevator Closure 

The control law for the overall pitch rate to elevator feedback is 

' {u-rfkj's + Mj <=»(=) (c-4) 

where 

Kq = KQAL 

Kq = KIQ 

Gg(s) = servoactuator transfer function 


Thus, as the effective feedback transfer function is 



Kq(s + 1/Te) ^ , 

(s + l/T„o) 


(C-5) 
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where 


-L = -1- + ^ 


The minus sign in Eq* C-5 is consistent with the assumption of negative 
feedback# For the BCS system as implemented on the PA-30 


— = 0.125 rad/sec 

^wo 


rp- = 3.70 rad/sec 

Kq 


Te 




3.83 = ~ rad/sec 

Kq 


Gs(s) 


1 

18s + 1 


The closure of the pitch rate loop is shown in the system surveys of 
Fig# C-1 for 80 and 110 kts. The Bode root loci (4) shown as dots in 
Figs. C-la and C-lb indicate the variation of the closed loop short 
period root with gain. Likewise the a root loci labeled (1) through (3) 
in Figs. C-la and C-lb show the variation of the closed loop real roots 
with gain. The roots for a given loop gain are indicated by the inter- 
section of the appropriate 0 dB line with the Bode loci. The roots cor- 
responding to the gain =* 0.25 rad/sec, are shown in Fig. C-1. The 
Bode loci in Figs. C-la and C-lb may be correlated with the conventional 
root loci shown in the generic sketch of Fig. C-1. 
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i.oa w(rad/sec) in.aa 
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b) Bode Roof Locus, ilOkts c) Generic Sketch Of Conventional Roof Locus 




The nature of the closed loop pitch rate dynamics is determined by 
the gain or, equivalently, the crossover frequency which is the 
frequency at which the 0 dB line for the chosen gain intercepts the 
Bode asymptote* It may be seen that, for the given values of the param- 
eters and i/Tg, the closure can be made in a K/s region which 
implies adequate phase margin and adequate low frequency gain margins* 
Thus, even though the short period and servo modes are highly coupled 
(especially at 110 kts) the closed loop characteristics should be 
acceptable - unless there is additional significant loss in phase margin 
due to time delays resulting from use of the ground computer. The char- 
acteristic equation with the pitch rate loop closed (and with implicit 
closure of the auto throttle) is 


A' = 


q^6e 

U‘^<S'p 



Kqd/Tg) 


1/Ta 

(1/Ta) 


M5^s2[1/Tq2) 


2 ^ ^ ^ ) [ ^ > ^ ] 


(C-6) 


It is of interest to consider the nature of the pitch rate response 
to inputs at the last summer before the servo. Inputs at this point may 
be considered to be pitch rate commands as denoted by ^q^^. Fig* A-1. 
Thus, the pitch rate feedbacks form a ^ pitch rate command system' with 
respect to q^ with the closed loop transfer function 
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qc 


1 


q>6j 


= 

qc 




Tas + 




A' 


TTTtTT ^ 

U/T^oJU/TaJIc'.'-tiJ 


(C-7) 


It may be seen from Figs. C-la and C-lb that the pole resulting from the 
washout time constant will essentially be driven into the attitude 
numerator l/T^^ at the nominal pitch rate closure. Thus 

/ -Ms (l /Two] 1 /Ta 

t ■ u',<j ‘xiTTn 


indicating that the pitch rate response to pitch rate command will be 
analogous to the open loop pitch rate response to elevator-command 
except that the attitude zero in Eq. C -8 is set by the washout time con- 
stant instead of the conventional value 1 /T 02 * 

4. a2 4e ^^SURE 


From Fig. A-1 , the a^ to elevator feedback transfer function is 
(negative feedback) 



(C-9) 


where 


« KINZ(37.2/57.3) 


The location of the accelerometer on the PA-30 was not specified in the 
available documentation, however, the analysis in Ref. C-2 is consistent 
with the fact that the accelerometer is near the instantaneous center of 
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rotation for elevator Inputs (ICR). Under this assumption and the fur- 
ther assumption that the elevator effectiveness in lift, may be 

neglected (so that the ICR is near the c.g.), then the accelerometer 
output response to elevator inputs for the speed constrained condition 
is given by Eq. C-10. 




u 

6 ^ 


nI 


N6x 


UoZwX6xM6eS _ ^oZyM^e 
XgjSAsp Asp 


The effect of the accelerometer feedback on the characteristic equa- 
tion with the pitch rate-to-elevator and auto throttle loops closed may 
be seen from a system survey of 



Ka/oZwM6e(^/T WO )(i/Ta) 


as sketched in Fig. C-2. The values of the accelerometer gain 
derived from Ref. C-2 are 


= 0.0768 > 0.1919 deg/g 


= 4.16 • 10“^ 1.04 • 10“^ rad/ft/secZ 
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\ 


w (rad /sec) —*~ 


X 

Ta 



K/s3 


\ 


Figure C-2. Asymptote Sketch of the Bode Root Locus for the 
Closure of the Loop (q>6g Closed) for 

the Existing System 


As shown in Fig. C-2, this is a very low value of accelerometer gain and 
thus this loop will have a negligible effect on the aircraft dynamics. 
This fact was also noted in the review of the linear analysis of 
Ref. C— 2 given in App. A. Consequently, the effect of the accelerometer 
feedback will be neglected for in the following analysis. 

5. SINK RATE TO ELEVATOR CLOSURE 

With the pitch rate loop closed and the accelerometer feedback 
neglected, the open loop transfer function for the h loop is 


# 

h_ 

>^3 




i A ±_ 

s q qc 


q+5e 


(C-12) 


For the short period equations of motion, the h response to pitch rate 
is given by 
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h . UoY . Up 

q s0 “ sCTe^s + 1) 

vhere 


T 



(-Zw 


+ 





Thus, Eq. C-12 becomes 


(C-13) 


UpZw^6e( 1/^wo ) 

s2(i/T92)[?'.u)A] * tl/TAJ 


and a system survey of this closure is sketched in Fig. C~3. For the 
crossover region defined by the nominal h gain at = 110 kts 



137Kh 

s 


(C-15) 


for 


Kh = 0.275 deg/fps = 4.8(10~3) rad/fps 


The sink, rate loop closure is made in the K/s region between the washout 
inverse time constant, and the airframe attitude zero, 1 /T 92 « 
The asymptotic crossover frequency is thus 137 = 0.66 rad/sec 
which will be the approximate bandwidth of the sink rate loop closure. 
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Figure C-3. System Survey of the h>6e Closure 
(q^6g Closed) for the Existing System 


It should be noted for comparison with the proposed system to be discus- 
sed later that achieving a good sink, rate loop closure requires that 
1 /T^q « 1 /T 02 » a situation which will be different in the proposed 
system. When comparing this closure to that in App. A, it should be 
noted that the situation here is somewhat different because the signifi- 
cant attitude lags of the rate-of-climb sensor have been neglected and 
the speed constraint has effectively removed the phugoid mode. 

6. SERVO MODE STABILITY 

Reference C-2 indicated that problems with limit cycle oscillations 
in the elevator were encountered in the PA-30 flight test of the BCS 
system. The analysis in App. A also indicates the sensitivity of the 
relative stability of the servo mode to changes in system parameters. 
Thus, a brief analysis of the servo mode will be presented here to 
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Identify the primary contributors to low damping ratio. For this 
analysis, a more complete servo model than the first order model used 
previously is required and the second order model given in Ref. C-2 as 


Gg(s) 



2 


(C-16) 


26,62 

" [0.70, 26.6] 


will be used. Since the pitch rate loop is the 'tightest' loop, it 
would be expected to be the primary cause of low servo mode damping 
ratio. This assumption may be tested by considering the expression for 
the characteristic polynomial with the pitch rate loop closed 


A' = 


.2 a 


Kqd/Ts) Ws 


(1/Two) 




(1/T02)s 


(C-17) 


Equation C-17 is equivalent to the expression in Eq. C-6 except for the 
use of second order servo model. At high frequencies near the servo 
mode, Eq. C-17 may be approximated as 


A' 


s^ - 


Kqa)gM§^s3 

[ ^s » “s J 


0.25 • 26.62 . _ 40.1 
(0)[0.70, 26.6] J 


^ s2(19.60) [0.464,19.0] 

[0.7, 26.6] 
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for 


Kq = 0.25 rad/ (rad/sec) 


(C-18) 


The closed loop servo poles given in Eq. C-18 check almost exactly with 
those given in Fig. A-4 of App. A which also includes the effect of the 
sink rate feedback closure with the ROL sensor lags. Thus, it may be 
seen that the effect of the pitch rate closure is dominant on the servo 
mode and that it can be expected that other lower bandwidth outer loops 
will not affect servo mode stability. The primary concern regarding 
servo mode stability will thus be effective time delays due to transmis- 
sion of the q>6g signal the airplane to ground computer and back to the 
airplane. 
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APPENDIX D 


ANALYSIS OF PROPOSED SYSTEM 


!• Design Concept For Beam Follower 

The basic design concept for the beam follower portion of the pro- 
posed autoland system makes use of the synthetic ILS beam derived from 
tracking radar data and presently used only as a pilot flight director 
display for manual landing. The approach taken here is to design a beam 
following system using conventional ILS system practice as presented in 
Ref. D-1 and D-2. For purposes of system analysis in the following sec- 
tions, it will be assumed the present PA-30/BCS schedules for sinkrate 
and airspeed will be used with the further assumption that beam tracking 
will occur in the decel region from an altitude of 950 ft to 150 ft with 
airspeed steadily decreasing from 110 kts to 80 kts. Landing flare 
using the variable t flare concept will be assumed to take place at a 
constant airspeed of 80 kts beginning at 150 ft altitude. These assump- 
tions are made to make the analysis concrete and to allow comparison 
with the existing BCS autoflare system. However, it will be shown that 
the basic system design and system gains can be used with variations in 
the altitude and airspeed schedules and furthermore, that the airspeed 
and altitude schedules can be considered independent. The primary con- 
cern in this connection is the effect of variations in airspeed during 
beam tracking. While the deceleration should always be slow enough with 
respect to the vehicle dynamic time constants that quasi-constant coef- 
ficient dynamic models may be used for analyses, there is at least a 
potential problem in the large variations in the effective gain through 
large changes in the dynamic pressure, e.g.. 


Mg ^ 

e 80 kts 
^elllO kts 



0.53 
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However, it will be shown in the following analysis that one set of sys- 
tem gains will work satisfactorily over the 80 to 110 kt speed range. 

A basic beam follower system concept is shown in Fig* D-1 and is the 
same as the Ref. D-2 design except that the pseudo-integral of pitch 
rate (from the pitch rate gyro) is used as a surrogate for pitch atti- 
tude, consistent with the mechanization of the existing PA-30/BCS 
system. 

It should be noted that the pitch rate feedback structure shown in 
Fig. D-1 is identical in form to that used in the existing PA-30/BCS 
autoflare system (Fig. A-1). No explicit accelerometer feedback loop is 
shown in Fig. A-1 and this is essentially consistent with the existing 
PA-30/BCS system given the negligible effect of the loop (see 

Appendix A analysis). It would, however, probably be desirable to com- 
bine the accelerometer and barometric rate-of-climb signals in a comple- 
mentary filter to produce a baro-lnertial sinkrate signal for the h 
feedback loop in Fig. D— 1. The primary difference in the sinkrate loops 
between the existing PA-30/BCS autoland system and the proposed beam 
follower is then the removal of the forward loop integrator from the 
existing system. 

2. Basic Design Considerations for Beam Follower 

The structure of the beam follower system shown in Fig. D-1 follows 
from consideration of the basic guidance, control and regulation 
requirements for the system. These are stated in Ref. D-1 and D-2 and 
will be repeated here for convenience: 

1) To establish and maintain the aircraft on a spe- 
cified spatial pathway or beam, e.g. the instru- 
ment low approach system ILS glidepath. 

2) To reduce flight path errors to zero in a stable, 
well damped and rapidly responding manner. 

3) To establish an equilibrium flight condition. 

4) To limit the speed or angle of attack excursions 
of this established equilibrium flight condition. 
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The regulation requirements are similar, i*e., 

1) To maintain the established flight path in the 
presence of disturbances such as gusts, winds, 
and wind shears 

2) To provide a degree of short term attitude sta- 
bility in the presence of disturbances. 

The requirements stated above relate primarily to the low frequency 
response of the aircraft /autoland system, i.e., they pertain to the 
desired spatial trajectory during beam tracking. However, achieving 
these basic guidance requirements requires an acceptable design for the 
inner attitude loops. The synthesis of the Fig. D-1 structure from the 
above requirements is analyzed in detail in Ref. D-2 and will be 
reviewed briefly here for perspective in the selection of gains for the 
proposed system. The most fundamental guidance requirement is that the 
aircraft establish itself on the beam when the system is engaged. In 
other words, the deviation from the beam, d, must ultimately become zero 
when the system input is an initial condition ^ 0. It is demon- 

strated in Ref. D-2 that this implies the need for a proportional, i.e., 
term in the beam deviation feedback loop to establish static stabil- 
ity with respect to the beam. 

It is also desirable that the system prevent beam deviation stand of 
errors, with minimal Iqng term error. This consideration leads to a 
requirement for an integral term in the beam deviation feedback (or even 
higher order integration paths depending on the complexity of the devia- 
tion command). For example, a change in glldeslope angle could be gen- 
erated by input of a ramp beam deviation command, d^. In the steady 

state, beam deviation error, d^, should go to zero: however, the 

measured slnkrate H will not be zero. Thus, to achieve a zero pitch 
rate command in the steady state, a non-zero value of h^, is required 
which in turn requires the presence of a Integrator in the beam 


A ... • 

Total slnkrate H = + h is measured where H is the nominal (non- 

o • • • 

zero) slnkrate for ideal beam following and h = d is a perturbation 
about the nominal. 
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deviation loop. However, since the beam following system will be phased 
out during flare in the proposed autoland system, the requirements for 
command following are less critical for this application. 

As noted in Ref. D-2, regulation against steady U-gusts imply a need 
for a proportion beam deviation feedback and regulation against horizon— 
tal wlndshear implies a need for an integral term, i.e., The 

requirements for regulation against W—gusts are more critical in that 
the integral term is required to regulate against steady gusts and that 
even with an integral term there will be a steady state error due to any 
shears component of a W~gust. W— gust regulation can be improved by 
selecting T^q properly. Ideally this is done by setting ^/'^wo 
below the short period frequency so that there is effective attitude 
(q/s) feedback at the short period mode, which is then washed out at 
lower frequencies. This allows the aircraft to weathercock at lower 
fr 0 quencies and thus unload W—gust induced angle of attack increments 
which would tend to cause the aircraft to deviate from the beam. 

•The above conclusions are derived from analyses in Ref. D—2 which is 
based on a 3 DOF analysis with unconstrained speed dynamics. While 
will be some difference in the dynamical details for the speed 
constrained case, the same basic system elements will still be required. 
Furthermore, given the uncertainties regarding the auto throttle design, 
it is desirable to use a system which is Independent of speed constraint 
as is the Fig. D-1 system. The inner slnkrate feedback loop approxi- 
mates a d feedback and thus provides a damping term in the beam devia- 
tion (higher frequency) dynamics. The pitch feedback inner loops pro- 
vide short period damping as well as short term attitude stability 
through the pseudo pitch attitude (pseudo integral of q) term. Finally, 
the inner pitch rate feedback provides a properly calibrated pitch rate 
command system for the outer path loops independent of uncertainties in 
the airframe pitch dynamics. 

3. q ♦ Se Loop Design 

The basic problem in design of the pitch rate feedback loop is the 
selection of the two time constants T^^ and Tg where 
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(D-1) 


Te = Kq/Kg 

Furthermore, It is desired to select one set of gains which are adequate 
over the 80 to 110 kt speed range. Figure D-2 shows a comparison of the 
Bode asymptotes at 80 and 110 kts for the airframe |q/6el under the 
speed constrained assumption and it may be seen that while there is 
almost a factor of two difference in dynamic pressure, the net effect on 
the response is not large. 

The basic considerations in the selection of the time constant Tg 
(which indicates the relative mix of effective pitch attitude and pitch 
rate feedbacks) may be seen as sketched in Fig. D-3. An upper limit on 
the value of the inverse time constant 


l/Tg < a>ca 


(D-2) 


is set by a desire to have the open loop Bode magnitude asymptote be 
K/s in the region of crossover for good open loop phase margin and 
consequently good closed loop damping for the short period mode. A high 
value of 1/Tg, approaching the crossover frequency, is desired, however, 
to maximize the low frequency gain margin, GM. This conflicts with a 
desire to place 1/Tg well below the crossover region to maximize the 
phase lead contribution at crossover. The desire for maximizing phase 
lead is generally not as critical as the need for good low frequency 
gain margin-unless the high frequency lags from the servo, etc. are very 
large. Thus, a near optimal selection for 1/Tg would be at or near the 
asymptotic crossover frequency, Uca* The value of 1/Tg = 3.7 rad/sec 
that is used for the existing PA-30/BCS system is a reasonable choice 
for the 80 to 110 kt speed range. Examination of the Bode plots for the 
existing PA-30/BCS system shown in Figs. C-la and C-lb of Appendix C 
indicate that while 1/Tg could be somewhat higher, this selection 
provides adequate low frequency gain margin and further allows some 
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u) (rad/sec) 


Figure D-2. Effect of Airspeed on Airframe 
|q/6el Frequency Response 



low high 


log cu 

Figure D-3 . Selection of l/Tg 
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reduction in the gain if needed to account for increased phase lags 
from time delays, etc. 

As noted previously in the discussion of guidance requirements for 
the beam following system, a primary consideration in the selection of 
the washout time constant is a desire to allow the aircraft to weather- 
cock at frequencies below the short period to reduce beam deviation due 
to W-gusts (i.e., windproofing). This implies setting 1/T^^ near but 
below the short period frequency. A further consideration is the effect 
of the washout on low frequency gain margin. As indicated in Fig. D-4, 
if 1/T^^ washout lag is placed near the attitude numerator I/T 02 , as 
shown in Fig. D-4b, there will be no additional loss of low frequency 
gain margin beyond that implied by the selection of 1/Tg. However, as 
the washout is moved above I/T 02 towards the short period there Is addi- 
tional loss in low frequency gain margin (Fig. D-4c). By way of compar- 
ison the washout used for the existing BCS system is shown in Fig. D“4a. 
While this provides additional low frequency gain margin it would 
detract from the windproofing goal. Furthermore, the selection of the 
low value of the washout lag in the existing PA-30/BCS system is criti- 
cal to the design of the sinkrate loop as shown previously in Fig. C-3. 
It will be shown shortly, this requirement is not necessary for the pro- 
posed system. Thus, consistent with the Ref. D-2 design philosophy, the 
washout inverse time constant is set as 


1 

^wo 


1 



80 kts 


1.48 rad/ sec < 


T 0^ 


IllO kts 


(D-3) 


Having set the two time constants, it only remains to select the 
loop gain K^. The primary consideration here is to use a high gain 
(i.e., high crossover frequency) to maximize the bandwidth of the inner 
loop as a pitch rate command system — since the inner loop bandwidth 
ultimately constrains the achievable outer loop bandwidth. The primary 
constraint on the pitch rate gain is destabilization of the servo mode 
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as noted previously in Appendix C. Since the phase margin in the cross- 
over region is relatively insensitive to the washout time constant, the 
situation is very similar to the existing pitch rate system. Thus, a 
of 0.25 rad/sec corresponding to the best gain noted in Ref. D-3 will be 
used for the proposed system. 

4. Closed Loop q/q^ Response 

With the effective inner loop dynamics viewed as a pitch rate com- 
mand system, the pitch rate response to pitch rate command is 


f 

qc 


A 




. 1 . 

(l/Two) ''Tas + 


1/Ta 


(1/Ta) 




(1/T^o)(1/TA)[c',(uA] 

(1/Two)(1/Ta) 

-M6^(1/Tq2) 

1/Ta 

U' »i%]- 

(1/Ta) 

40.1(1.98) 

18 

[0.657, 14. IJ 

(4.45) 

21.8(1.48) 

18 

[0.889, 8.99] 

(6.01) 


110 kts 


80 kts 


(D-4) 
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5* Relation of the Inner Attitude and Outer Path Loops 


Using the speed constrained sinkrate-to-pitch rate response transfer 
function, Eq. C-13, Appendix C, the sinkrate-to-pitch rate command 
transfer function is 


i ^ M 1 ^o 

q s 9 s(l02S + 1) 

i- i • a-(s) 

qc s(T 02S + 1) qc 

Uq/'^02 + 1/T02) 1/T^ 

s(s + 1 /T 02 J WtIT 


Noting that I/T 92 == 

^ " TnTrxrrrr^ 


From Eq, D-5 it may be seen that the sinkrate-to-pitch rate command 
response is 'K/s' out to the closed loop actuator break point 


1/Ta 


4.45 rad/sec at 110 kts 
6.01 rad/sec at 80 kts 


(D-6) 
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which is above i/T 02 < 1*9 rad/sec. Because of the strong actuator- 
short period coupling, the sinkrate loop bandwidth might even be 
improved by reducing although this should not be necessary. When the 
succeeding beam deviation loops, d and Jd to elevator are closed, there 
will thus be a hierarchy of K/s closures. 

6* Selection of Path Loop Gains 

The control law implied by Fig. D-1 is 


<5e 


Kd + KdS 
s(TfS + 1) 


d + K^d 


+ Gq®q 


Gg(s) 


Kd[l/Tdi)(l/Td2)(l/Td3) 

^ITTtTI 


d + GqSq 


G„(s) 


(D-7) 


As noted in Ref. D-2, the beam deviation inverse time constants are 
approximately related to the system gains as 


1 ^ 51.' _L_ ^ Ik JL 1 J_ 

Tdi ~ *^d ’ Td2 Kj * Td3 “ Tf 


(D- 8 ) 


Selection of the beam filter characteristics which define the time con- 
stant, depend on the nature of the beam noise which is unknown for 
the synthetic ILS beam. However, it may be assumed that the time con- 
stant cancel 1/Tp. Thus, only the tin^ constants 1/T^ and 

1/Tj need be selected. As noted in Ref. D-2 these time constants are 
^2 

ordinarily related to the phugoid but under the assumption of speed con- 
straint, we merely need to select the inverse time constants to be well 
below the beam deviation loop crossover frequency to achieve a wide K/s 
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region with near maximum phase lead from these terms. Therefore, fol- 
lowing Ref. D— 2 we select 


0.3 rad/ sec 

(D-9) 

0.15 rad/ sec 

The beam follower loop closures are summarized in the Bode root loci 
of Figs. D-5a and D-5b at 110 and 80 kts. For comparison, a generic 
sketch of the conventional root locus for the beam follower is shown in 
Fig. D-5a. It may be seen that for any speed in the 80 to 110 kt range, 
there is a broad K/s region between 1/T^j^ and the actuator time con- 
stant 1 /Ta with sufficient phase margin to close the beam follower loop 
at crossover frequencies approaching 2 rad/sec. On this basis, a loop 
gain, Kd = -42 dB or 0.0080 rad/ft/sec, should provide good beam follow- 
ing performance. It should be noted that there is a conditional low 
frequency instability in this system which implies a lower limit on the 
gain Kd well as an upper limit. These two limits correspond roughly 
with the extent of the K/s region and do not compromise the design. The 
low frequency conditional instability, which would occur only if the 
gain Kd were too low, is a consequence of the speed constraint and does 
not occur for the unconstrained 3 DOF analysis in Ref. D-3. 

7. Variable x Flare System Design 

A generic mechanization of the variable r flare law is shown in 
Fig. D— 6 . There are several possible ways to obtain a measure of ground 
speed Vq in wind which will affect details of the actual mechanization, 
but will not affect selection of the basic gains shown in Fig. D— 6. 
Figure D-6 produces a variable x flare law consistent with Eq. 24 which 
(in a rearranged form) is 


1 . Kd 

^d2 Kd 


1 . M 

Tdi Kd 
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<j(rod/sec) 


Conditional 

Instability 







Figure D-6. Generic Variable r System Mechanization 


dh , Vg 

dt tqVgq 


(h + hg) = 0 


(D-10) 


The above ordinary differential equation is the ideal (hg = 0) for sum- 
mer A in Fig. D-6. That is, the signal through the upper forward path 
is proportional to altitude, h, where 


^ 1 * ■ 


(D-11) 


and the lower forward path is proportional to sinkrate dh/dt. The mech- 

# 

anization shown in Fig. D-6 causes the baro-inertial sinkrate, \iQ-p used 
in beam following to be washed out in flare with time constant Tp and 
replaced by Che derivative of the radar altimeter signal, i.e., washed 
out hr, which is "runway referenced." This system can thus be considered 
to be runway referenced within 3 Tp seconds after flare initiation. 
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Thus, the primary consideration in the selection of the washout time 
constant, Tp is that it be much shorter than the flare law time con- 
stant T* 

The control law derived from Fig. D-6 is 


Kh 


KhVGhR + (tfs% i)^R (t'fs + l) + hg 


(D-12) 


If it is assumed for basic analysis that the radar altimeter and baro- 
inertial rate of climb sensor characteristics are ideal (unity gains) 
then 


= Kh 


KhVc + 


s^Tp 


Tps +1 Tps + 1 


h + Khhg 


= Kh 


s2 + (K^Vg + 1/Tp)s + K^Vg/Tf 
( s + 1/Tp) 


h + K,;hB 


(D-13) 


Since Kj^Vq = 1/x 


Kh(s + l/x)(s + 1/Tp) 
(s + 1/Tp) 


h + KhhB 


= + l/x)h + Kj^hB 


(D-14) 


Thus, from Eq. D— 14 it may be seen that the Fig. D— 6 mechanization pro- 
duces the variable x control law as required. 
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The same inner loop pitch rate command system used for the beam fol- 
lower may be used for the flare law and thus, the open loop sinkrate-to- 
sinkrate error transfer function is (from Eq. D-5) 


^ • K 

qc “ s(1/Ta)[c',<] 


(D-15) 


where K = 


Equation D-14 implies a feedback transfer function 
Gh® = -Kh(s + 1 /t) 


(D-16) 


The open loop transfer function for the flare law is thus 


hi 

qc 


OL 


- KhK(s + 1/x) 
s qc ~ s2(l/TX)[c'.a»'] 


(D-17) 


which is shown in the bode root locus plot of Fig* D-7 for an 80 kt 
flare with the nominal (zero wind) t = 19*8 sec corresponding to the 
existing PA--30/BCS law* It may be seen from Fig. D-7 that there is a 
broad K/s region with potential for a 2 to 3 rad/sec crossover for the 
flare law. Unlike the beam follower no conditional instability occurs 
in this case because there is no integral feedback term. It may be seen 
that the flare inverse time constant is much lower than the desired 
crossover frequency and thus the crossover region is quite insensitive 
to 1/t over the ±30 kt steady wind range where 

0.032 ^ < 0.069 (D-18) 
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I Closed loop root 
' |G (jw)| and 4 G( ju) open loop 

|G(jw)|ond 4G(jcj)open loop, asymptote 

|G(-cr)|left half plane (stable) root locus 

• • • • Complex plane root locus 



Figure D-7 . Closure of Flare Low Loop in Figure D-6 
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Thus, the variation of x over the nominal windspeed range should not 
effect the flare law loop closure, A selection of = -38 dB = 0.0126 

rad/ft/sec will give an asymptotic crossover frequency of about 2 rad/ 
sec. 
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APPENDIX E 


CIRCULAR ORBIT ALGORITHM 


A. TIME HISTORIES 

Figures E-1 to E-5 present time histories for several values of air 
speed and wind speed (dimensions in ft, sec, rad). Variables are XI 
X, X2 = y, X3 = \|;. 
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Figure E-2. Circular Algorithm, U = 400, = 100 
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Figure E-3. Circular Algorithm, U = 300, = 100 
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Figure E-4. Circular Algorithm, U = 250, = 100 
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Figure E-5. Circular Algorithm, U = 200, = 100 
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APPENDIX F 


DEAD RECKONED WING POINTING ALGORITHM 
A. TIME HISTORIES 

Figures F-1 to F-3 present time histories for several values of air 
speed and wind speed (dimensions in ft, sec, rad). Variables are XI 
X, X2 = y, X3 = ip. 
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Figure F-1 . Wing Pointing Algorithm, U = 500, V^.= 0 
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Figure F-2. Wing Pointing Algorithm, U = 500, = 100 
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Figure F-3. Wing Pointing Algorithm, U = 250, = 100 
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